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Abstract 



We report the observation of bottom-charmed mesons B c in 1.8 TeV pp 
collisions using the CDF detector at the Fermilab Tevatron. The B c mesons 
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were found through their semileptonic decays, Bf- — > J/ip£^X. A fit to the 
J/ip£ mass distribution yielded 20.4^'g events from B c mesons. A test of 
the null hypothesis, i.e. an attempt to fit the data with background alone, 
was rejected at the level of 4.8 standard deviations. By studying the quality 
of the fit as a function of the assumed B c mass, we determined M(B C ) = 
6.40 ± 0.39 (stat.) ± 0.13 (syst.) GeV/c 2 . From the distribution of trilepton 
intersection points in the plane transverse to the beam direction we measured 
the B c lifetime to be t(B c ) = 0.46±g;i| (stat.) ± 0.03 (syst.) ps. We also 
measured the ratio of production cross section times branching fraction for 
B+ -> J/^t+v relative to that for B + -> J/^K+ to be: 

^(b)-Ir^j/^k) = °" 132 ^^ (stat) ± °- 031 (syst - } (lifetime) 

PACS numbers: 14.40.Nd, 13.20.He, 13.30.Ce, 13.60.Le, 13.87.Fh 
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I. INTRODUCTION 



The meson is the lowest-mass bound state of a charm quark and a bottom anti- 
quark.0 It is the pseudoscalar ground state of the third family of quarkonium states. Since 
the B c has non-zero flavor, it has no strong or electromagnetic decay channels, and it is the 
last such meson predicted by the Standard Model. Its weak decay is expected to yield a large 
branching fraction to final states containing a J/ip a useful experimental signature. 

Non-relativistic potential models are appropriate for the B c , and they predict its mass. 
Kwong and Rosner ||] estimate M(B C ) to be in the range 6.194-6.292 GeV/c 2 . Eichten and 
Quigg || discuss four potentials that yield values in the range 6.248-6.266 GeV/c 2 . In these 
models, the c and b are tightly bound in a very compact system. These authors describe a 
rich spectroscopy of excited states, which make this the "hydrogen atom" or, perhaps, "the 
mu-mesic atom" of QCD. 

We expect the full decay width of the B c to consist of three major contributions, T = 
T b + T c + T bc , which are, respectively, 

• b — > cW + with the c as a spectator, leading to final states like (J/ipn), (J/i/jiu); 

• c — > sW + , with the b as spectator, leading to final states like (B s n), (B s lv); 

• cb —> W + , annihilation leading to final states like (D* K), (t u T ) or multiple pions. 

Since these processes lead to different final states, their amplitudes do not interfere. In 
the simplest view, the c and b are free, so annihilation is suppressed, and the total width 
is just the sum of the c and b total widths, with c-decay dominating. Approximating this 
by F(B C ) = T(D°) + T(B°) yields t(B c ) « 0.3 ps |§. When annihilation, phase space 
considerations (which reduce the relative importance of the c contribution) and other effects 
are included, the predictions increase to the range 0.4-0.9 ps [p],|T|-|lCi|j . Quigg [T1J emphasizes 



the relatively large ratio of the binding energy to charm-quark mass and the effect on r& c 
of the compact size of the cb system, where the pseudo-scalar decay constant is expected 
to be fs c ~ 500 MeV. He predicts lifetimes in the range 1.1—1.4 ps, with Tb as the largest 
contribution. Thus, a B c lifetime measurement is a test of the different assumptions made 
in the various calculations. Several authors have also calculated the B c partial decay rates 
to semileptonic final states |]-fll|l2| . 

In perturbative QCD calculations of B c production using the fragmentation approx- 
imation, the dominant process is that in which a b is produced by gluon fusion in the 
hard collision and fragmentation provides the c [p~3| p~Tf| . A full a\ calculation shows that 



fragmentation dominates only for transverse momenta large compared to the B c mass, i.e. 
Pt 3> Mb c c jrj3|. This calculation provides inclusive production cross sections along with 



distributions in transverse momentum and other kinematic variables. 

There have been several experimental searches for the B c meson. In e + e~ collisions at the 
Z resonance at LEP, 90% confidence level (C.L.) upper limits have been placed on various 
branching- fraction products by the DELPHI collaboration []18], the OPAL collaboration [19 



1 References to a specific state imply the charge-conjugate state as well. 
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and the ALEPH collaboration ||20|| . In Sec. |VIII| , we compare these limits with our result. 
OPAL reported one event in the semileptonic channel where the background was estimated 
to be (0.82 ± 0.19) event, along with two Bf — > J/ipir^ 1 candidates with an estimated 
background of (0.63±0.20) events. The mean mass of the latter two candidates is (6.32±0.06) 
GeV/c 2 . ALEPH [^0| reported one candidate for — > J/ijj^v^ with a low background 
probability and a J/ip iM mass too high to be explained by a light B meson. A prior CDF 
search placed a limit on the production and decay of the B c to J/ip and a charged pion [ET 



We report here the observation of B c mesons produced in a 110 pb -1 sample of 1.8 TeV 
pp collisions at the Fermilab Tevatron collider using the CDF detector. We searched for the 
decay channels B c — > J/ip \iX and B c — > J/tpeX with the J / ip decaying to muon pairs.0 
Even the lowest prediction for the B c lifetime @ implies that a significant fraction of J/ip 
daughters from B c would have decay points (secondary vertices) displaced from the beam 
centroid (primary vertex) by detectable amounts. The existence of an additional identified 
lepton track that passes through the same displaced vertex completes the signature for a 
candidate event. We have identified 37 events with J/ip£ mass between 3.35 GeV/c 2 and 
11.0 GeV/c 2 . Of these, 31 events lie in a signal region 4.0 GeV/c 2 < M{J/ip£) < 6.0 
GeV/c 2 . 

The most crucial and demanding step in the analysis is understanding the backgrounds 
that can populate the mass distribution [^2]f23|. We attribute any excess over expected 



background to production of the B c , the only particle yielding a displaced-vertex, three- 
lepton final state with a mass in this region. The bulk of the background arises from real 
J/ijj mesons accompanied by hadrons that erroneously satisfy our selection criteria for an 
electron or a muon or by leptons that have tracks accidentally passing through the displaced 
J/ip vertex. 

In the sections that follow, we begin with a very brief discussion in Sec. [TJ of some parts 
of the CDF detector, particle identification, and identificaton of J/ip through its decay to a 
muon pair. Following this, we describe our selection criteria for tri-lepton events (Sec. JTJ). 
our calculation of the number of background events in the signal region (Sec. |TV|) , and the 
validation procedures to establish the accuracy of that calculation (App. [B]). 

Section [V] describes the procedures we used to establish the existence of the B c contri- 
bution to our sample of candidates. The background calculations and the mass distribution 
of the J/ipl data sample were subjected to a statistical analysis from which we calculated 
the B c contribution to the signal region. We describe first a simple "counting experiment" 
calculation for events in this region. However, we base our claim for the existence of the B c 
on a likelihood fit that exploits information about the shape of the signal and background 
distributions in the mass range 3.35-11.0 GeV/c 2 , which we call the fitting region. The B c 
contribution to these data is 20.4^g' 2 events. The null hypothesis is rejected at a level of 4.8 
standard deviations, i.e. the probability that the background could fluctuate high enough 
to explain this excess is less than 0.63 x 10~ 6 . 

In Sec. |VlL by studying the quality of the fit as we varied the assumed B c mass, we 



J/ipiv 



2 Because of the large partial widths for B c 
dominate B c — > J/ip£X , and we often refer to them simply as B, 
Sections VII and |VII| we discuss this further. 



12], we assume that these modes 
> J/ip£v or J/ip£. In 
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obtained an estimate of M(B C ). In Sec. VII], we describe our measurement of the B c lifetime, 



and in Sec. [VIII| we describe our measurement of the cross-section times branching-fraction 
ratio: 

a(B+) ■ Br(B+ -+ J/j)l + v) 
a(B+).Br(B+^J/^K+) 

We chose this form because many of the uncertainties cancel in the ratio. 
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II. DETECTOR AND PARTICLE IDENTIFICATION 



We collected the data used in this analysis at the Fermilab Tevatron Collider with the 
Collider Detector at Fermilab (CDF) during the 1992-1995 run. The integrated luminosity 
was 110 pb -1 of pp collisions at y/s = 1.8 TeV. We have described the CDF detector in detail 
elsewhere |]2^ , |25[| . We describe only those components that are important for this report. 



The events we sought, B c — > J/ifj£v where J ftp — > have a very simple topology: 

three charged particle tracks emerging from a decay point displaced from the primary inter- 
action point. For each track, the momentum must be known, along with its identity, fi or e. 
Below we describe the charged-particle tracking system, the electron identification system, 
the muon identification system, the real-time triggers, and J/ip identification. 



A. Charged Particles 

Our cylindrical coordinate system defines the z axis to be the proton beam direction, 
with cf) as the azimuthal angle and r as the transverse distance. Three tracking subsystems 
detect charged particles as they pass through a 1.4 T solenoidal magnetic field. We discuss 
them in order of increasing distance from the beam axis. 

• The silicon vertex detector, SVX, provides r— <p information with good resolution close 
to the interaction vertex. It consists of four approximately cylindrical layers of silicon 
strip detectors outside the beam vacuum pipe and concentric with the beam line. The 
active area of silicon is centered within the overall CDF detector and extends 25.5 cm 
in each direction along the beam line. The four layers of detectors are at radii of 3.0, 
4.2, 5.7, and 7.9 cm |^6|,^7|. The strips are arranged axially, and have a pitch of 60 



/im for the three innermost layers and a pitch of 55 /im for the outermost layer. 

• A set of time projection chambers provided r—z information that was used to determine 
the event vertex position in z, which serves as a seed in the reconstruction of tracks 
in the r-z view in the drift chamber described next. 

• The central tracking chamber (CTC) is an 84-layer cylindrical drift chamber, which 
covers the pseudorapidity interval \r]\ < 1 (where rj = — ln[tan(#/2)] and 9 is the polar 
angle with respect to the proton beam direction). It consists of five superlayers of axial 
sense wires interleaved with four small-angle stereo superlayers at an angle of about 3° 
with respect to the axial wires. In each axial (stereo) superlayer there are twelve (six) 
cylindrical layers of sense wires. The efficiency for track reconstruction is about 95% 
and independent of pt for tracks with p? > 0.5 GeV/c. From the reconstructed tracks, 
we used charge deposition from hits in the outer 54 layers of the CTC to measure the 
specific ionization (dE/dx) of particles with about 10% uncertainty. This enabled us 
to determine the relative ir/K/p contributions in background calculations. Specific 
ionization was also used as one of the electron identification criteria. 

The combined data from SVX and CTC, required for all tracks in this analysis, have a 
momentum resolution Spx/pr = [(0.0009 x p T ) 2 + (0.0066) 2 ] 1//2 , where pt is in units of 
GeV/c, and the average track impact parameter resolution is (13 + (40/pt)) H m relative to 
the origin of the coordinate system in the plane transverse to the beam |f26f. 
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B. Electron Identification 



Electrons were identifed by the association of a charged-particle track with p? > 2 GeV/c 



and an electromagnetic shower in the calorimeter |25 |. The central (\rj\ < 1.1) calorimeter 
is divided into towers that subtend 15° in azimuth and 0.11 units of pseudorapidity. Each 
tower has two depth segments, a nineteen-radiation-length electromagnetic compartment 
(CEM) and a hadronic compartment. 

The track must project sufficiently far from a tower boundary that the energy deposition 
by an electromagnetic shower would be largely contained within a single tower. The energy 
E observed in the CEM tower must be roughly consistent with the momentum p of the 
track, viz., 0.7 < E/pc < 1.5, and we require that the energy in the hadron compartment of 
this tower be less than 10% of that found in the CEM. 

Information from other detectors further improves electron identification. The value 
of dE/dx measured in the CTC must be consistent with that expected for an electron. 
Pre-radiator chambers located between the magnet coil (one radiation length thick) and the 
CEM must show a signal equivalent to at least four minimum-ionizing particles. Proportional 
chambers with both wire and cathode-strip readout are located in the CEM at a depth of 
six radiation lengths. The shower profile observed in orthogonal views in these chambers 
must be consistent in pulse height, shape, and position with those found for electrons. 

Real electrons can arise from photon conversions to e + e~ pairs, including internal con- 
versions in 7T° — > 7e + e~. These can be identified and rejected when the candidate electron, 
paired with an oppositely charged track in the event, is kinematically consistent with the 
hypothesis 7 — > e + e~. However, such tracks were useful in direct measurements of our 
electron identification efficiency. 



C. Muon Identification 

Muons from J/ip decay were identifed by matching a charged-particle track with 
Pt > 2GeV/c to a track segment found in the muon drift chambers that lie outside the 
central calorimeter. The calorimeter presents five interaction lengths for \r]\ < 0.6 (CMU 
detector) and six to nine interaction lengths for 0.6 < |?7| < 1.0 (CMX detector). Within the 
uncertainty introduced by multiple Coulomb scattering, we required the charged-particle 
tracks found in the CTC and SVX to project to the track segments in these drift chambers 
within three standard deviations. 

We refer to the muon produced directly in the B c semileptonic decay as the "third muon," 
and we apply stricter requirements to identify it | 25f| . The transverse momentum of the third 



muon was required to exceed 3 GeV/c. A third muon must project to a track segment in the 
CMU, and for further suppression of backgrounds must pass through an additional three 
interaction lengths of steel to produce a track segment in a second set of drift chambers 
(CMP detector). These chambers cover about two-thirds of the solid angle for \r]\ < 0.6. 
Above 3 GeV/c, the efficiency for a muon track to match track segments in both the CMU 
and the CMP is independent of pr- 
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D. J ftp Selection 



The CDF detector includes a three-level, real-time trigger system with options that 
can be used to select events appropriate for a wide range of physics topics. In order to 
ensure consistent treatment for — > J/tp fi + X decays, — ► J/ipe + X decays, and the 
B + — > J/ip K + decays used for the cross-section normalization, we required that the muons 
from the J/ip decay satisfy the di-muon trigger selection requirements. 

The Level-1 trigger identified muon-chamber candidates by requiring a coincidence be- 
tween two radially aligned muon chambers. Our di-muon trigger required two such coinci- 
dences. 

The Level-2 di-muon trigger combined the muon candidates with information from a fast 
track processor that identified tracks from CTC data |28||. For the first 19.4 pb _1 of data 



collected, we required a single match between a muon chamber coincidence and a CTC track 
with pt > 3 GeV/c. The upgraded trigger system used for the remaining data required two 
such matches for tracks with p? > 2 GeV/c. Curves of the pt thresholds for the fast track 
processor and for the muon chambers can be found in Ref. [ 29f| . 

The Level-3 di-muon trigger was a preliminary event reconstruction in which we required 
charged muon candidate pairs with a mass, determined from CTC information only, between 
2.8 and 3.4GeV/c 2 . 

Subsequent offline processing performed a comprehensive search for all muon candidates 
in the event. For consistent treatment of the several decay modes described above, we 
required that the muons used to search for J/ip candidates were identical to those that 
triggered the event. We also required that both muons pass through the SVX. 

We performed a \ 2 fit to the track parameters for pairs of oppositely charged muons 
subject to the constraint that they had a common origin ||29|| . The di-muon mass was 
unconstrained. We required the \ 2 probability of the fit to exceed 1%. The resulting di- 
muon mass distribution is shown in Fig. [I]. The mean mass resolution is 16 MeV/c 2 . We 
required di-muon candidates selected by the offline programs for the B c analysis to be within 
50 MeV/c 2 of the world average J/ip mass of 3096.9 MeV/c 2 |0 . 
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III. EVENT SELECTION 



To identify B c candidates, we searched for events with a third track that originated at the 
J/ip decay point. We subjected the three tracks to a \ 2 fit that constrained the two muons 
to the J/ip mass and that constrained all three tracks to orginate from a common point. We 
accepted events for which the fit probability satisfied P(x 2 ) > 1%- To the resulting samples 
of J/ip + track, we applied further geometric and particle-identification criteria for selecting 
J/ip e and J/ip n events and a kinematic test for selecting J/ip K events.^ 

The third track for most events was a pion or a kaon.[] The fitting program corrected 
individual tracks for ionization losses. Consequently, the fit results had some slight sensitivity 
to the mass assumed for the third track. For studies aimed at identifying events with a 
specific third particle (e*, ^ or K^) we used the appropriate mass. For generic J/ip + 
track studies we used the muon mass. 



A. J/ip + track Decay Vertex Position 

The di-muon fit described in Sec. [Tl] constrained the daughter tracks from J/ip — > n + fi~ 
to come from a common point in space based on information from the CTC and SVX. When 
fitting the two muon tracks of the J/ip and the additional track, we required all three tracks 
to come from the same vertex. However, the high-resolution information from the SVX 
provides no longitudinal (z) coordinate. Thus, we measured the displacement between the 
beam centroid and the J/ip decay point in the transverse plane. The uncertainty in the 
displacement is typically about 55 /mi, and the uncertainty in the position of pp collision 
which produced the J/ip is 23 /im [p3 . 



L xy is the distance between the beam centroid and the decay point of a B c candidate 
projected onto a plane perpendicular to the beam direction and projected along the direction 
of the B c in that plane. A measure of the time between production and decay of a B c 
candidate is the quantity ct*, defined as 

M(J/^i).L xy (J/^£) 

IptOW)I [) 

where M(J/ip£) is the mass of the tri-lepton system and pt{J/iP£) is its momentum trans- 
verse to the beam. The average uncertainty in the measurement of ct* is approximately 25 
/im. In order to reduce backgrounds involving prompt J/ip production, we required ct* > 
60 /im for all candidates in the analysis of the B c signal significance. For the subsequent 



lifetime analysis (Sec. IVTJ), this requirement was modified. 



3 Differences in the criteria for identifying muons and electrons yielded different acceptances and 
backgrounds for the two decay channels. However, wherever it was possible to adopt common 
procedures for the two channels, we did so. 

Preliminary studies of dE/dx for the this sample of tracks showed the contribution from protons 
and antiprotons to be negligible and it was assumed to be zero thereafter. 
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B. J/ip K Identification 



The B — > J/ip K final state has no undetectable particles and can be reconstructed fully 
to calculate the mass of the parent B meson. We determined the mass for each J/ip + track 
combination under the hypothesis that the track corresponded to a kaon. 

Figure |2| shows the J/ipK mass distribution. The results from this particular data 
sample were used to normalize the measurement of the product of the B c production cross 
section and the B c — > J/ipiv branching fraction described in Sec. |V11J| . Events for which 
M(J/ip K) was within 50 MeV/c 2 of M{B) = 5.2789 GeV/c 2 were designated as B -> J/ip K 
and removed from the sample of candidates for B c — > J/ip I v. With different sets of selection 
criteria, the J/tp K sample was used to check the calculation of the probability for a kaon 
to be falsely identified as a muon (Sec. |1V A 1[ ) and to normalize Monte Carlo calculations 
of backgrounds from BB pairs (Sec. |IVD 



C. J/ip + Lepton Identification 

Figures 0(a) and 0(a) are histograms of the J/ip + track mass for combinations that 
passed the requirement P(x 2 ) > 1% described above. We required third tracks to have an 
opening angle less than 90° relative to the J ftp direction. This reduced the amount of BB 



background discussed in Sec. |IV|. The B — ► J/ip K events excluded from the J/ip + track 
sample populate a very narrow region of M( J/ip+track) in Figs. and 0. 

In the likelihood analysis described in Sec. [V|, the widths of the mass bins are not uniform. 
In Fig. |H and in subsequent figures containing mass histograms the bin boundaries are 
indicated by tick marks at the top of each figure. Most bins are 0.3 GeV/c 2 wide. We 
confined the effects of the excluded events near M(J/ip K) to one 0.15 GeV/c 2 bin, which 
is clearly visible in the figures. We also adopted wider bins at high masses where the 
event population is low. We chose the vertical scale so that the number of events per 0.3 
GeV/c 2 is equal to the number of events per bin for most bins. This makes explicit the 
statistical significance for the candidate distributions in Figs. 0(b) and f|(b). The event 
count is displayed for the two bins in Figs. 0(b) that had to be scaled. 

With an assumed B c mass of 6.27 GeV/c 2 , Monte Carlo simulations (App. [A]) reveal 
that 93% of the the tri- lepton masses reconstructed for J/ip ji and J/ip e decays will fall in 
the range 4.0 to 6.0 GeV/c 2 . We refer to this as the signal region. When we apply the 
muon identification criteria to events in Fig. 0(a), we obtain the mass distribution shown in 
Fig. f|(b), in which 12 of the 14 events lie in the signal region. When we apply the electron 
identification criteria described earlier to events in Fig. 0(a), we obtain the mass distribution 
shown in Fig. 0(b), in which 19 of the 23 events lie in the signal region. 

The distributions shown in Figs. 0(a) and 0(a) have many events in common because 
most with tracks that satisfy the muon and geometric criteria also have tracks that satisfy 
the electron p? and geometric criteria. Figures 0(b) and 0(b) have no events in common. 

The two candidate mass distributions contain irreducible backgrounds from various 
sources over the entire mass range. There are 37 candidates, of which 31 lie in the signal re- 
gion. Our principal task was to understand the shape and normalization of the backgrounds 
over the whole range of masses. We then determined their contributions to the signal region 
and established the size and significance of a B c contribution to that region. 
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D. Efficiencies 



The analyses described in the following sections required the relative values for the follow- 
ing efficiencies: e e = e(B c — > J/ipeX), £ IX = e(B c — > J/ippiX) and e K = e(B ± — > J/ip K^). 
We used a Monte Carlo program (App. |A|) to study the response of our detector and re- 
construction programs to each of these processes. All Monte Carlo events were subjected to 
the same requirements as the data. Among these requirements we emphasize ct* > 60 jum 
and M(J/ip£) in the range 3.35 to 11.0 GeV/c 2 . In order to eliminate shared systematic 
uncertainties, such as those associated with J/ifj detection, triggering and reconstruction, 
we used only the ratios of these efficiencies: 

R e = £e = 0.58 ± 0.04 (2) 

£ 4- £ 

c e ~ c fj, 

RK = £e_ = Q 244 ± Q33 ^ 

e K 

The principal differences between the efficiencies for J/ip e and J/ip fi are the larger 
geometric acceptance for the electron identification relative to that for muon identification, 
electron isolation requirements in the calorimeter, and the different pr thresholds: 2.0 GeV/c 
for electrons and 3.0 GeV/c for muons. 

The uncertainties in e e and that do not cancel in R e come from differing particle 
identification procedures for electrons (10%) and muons (5%), uncertainty in the Monte 
Carlo calculation (10%), and model dependence (App. [A]) due to the differing p T thresholds 
for muons and electrons (5%). This model dependence arises from uncertainty in the pt 
spectrum for B c production. As a check of our B c production model, we show in Fig. |5| 
the tri-lepton px distribution for the 31 candidate events in the signal region compared to 



those for simulated B c events and for calculated backgrounds (Sec. 0). There are no major 
differences in shape among the three distributions. 

The uncertainties in R come from Monte Carlo statistics (4%), uncertainties in the 
model (App. |A|) for production p? spectra (5%) and in the fragmentation parameter (2.3%), 
uncertainties in the detector (5%) and trigger (4%) simulations, and uncertainty in the 
electron identification (10%). 

We calculated the efficiencies for B c decays assuming a B c lifetime cr = 120 /mi. Lifetime 
effects cancel in R e but not in R K . R K scales as the number of B c that survive the 60 /im 
threshold in ct*, i.e. 



R K (cr) = R«(U0»m) { 1 K)CT \ ( l ; 

P V (1/^)120 Atm 

where (1/K)ct is the effective mean decay length, and the average correction factor is 
(l/K) = 0.88 ± 0.02. (See Sec. ^T|.) 



60 (im 
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IV. BACKGROUND DETERMINATION 



Backgrounds in the sample of B c candidates can arise from misidentification of hadron 
tracks as leptons {i.e. false leptons), from random combinations of real leptons with J /if) 



mesons, and from incorrectly identified J/ip candidates f22| , |23 



We describe three sources of false lepton identification. 

• The third track is a kaon or pion that has passed through the muon detectors without 
being absorbed. We call this "punch-through background." 

• The third track is a kaon or pion that has decayed in flight into a muon in advance of 
entering the muon detectors. We call this "decay-in-flight background." 

• The third track is a kaon or pion that has been falsely identified as an electron. We 
call this "false electron background. "Q 

Random combinations arise from the following sources: 

• External or internal conversions, i.e. electrons from photon pair-production in the 
material around the beam line or from Dalitz decay of tt°. Electrons from these sources 
that escape identification as conversions are called "conversion background." 

• KB that has decayed into a J/tp and an associated B that has decayed semileptonically 
(or through semileptonic decays of its daughter hadrons) into a muon or an electron. 
The displaced J/ ip and the lepton can accidentally appear to originate from a common 
point. We call this U BB background." 

Table | (Table |I|) summarizes the results of the data and background for the muon and 
electron channels in the signal (fitting) region defined in Sec. |. The procedures used to 
obtain these results are described in the remainder of this section. We have also conducted 
studies to verify the accuracy of our background calculations, applying them to independent 
data samples where they can be checked against direct measurements. These studies are 
described in App. [B| 

A. False Muon Backgrounds 

1. Punch- Through Background 

One of the backgrounds that can mimic a B c — > J/ip fiX event results when a n ± or 
or one of the particles in the resulting shower is not completely contained in the calorimeter 



5 As stated in Sec. Ill, we made the conservative assumption that the hadron tracks are all from 
mesons. Protons do not decay in flight. They have an interaction cross section higher than that 
for mesons and, therefore, a lower punch-through probability. Abandoning this assumption would 
lower our estimate of false muon backgrounds by a fraction of an event. The assumption does not 
apply to our procedure for estimating false electron backgrounds, which was validated with jet 
data containing a mix of mesons and baryons (App. ^). 
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and CMP steel. This can cause the original track to be misidentified as a muon. Although 
the probability for this is about 1 in 500, a large number of events have tracks that meet 
the fiducial requirements, which offsets the low punch-through probability. Such tracks can 
be reconstructed with a J/ip to mimic a B c decay. 

We used a model of the distribution of material in the CDF detector and the absorption 
cross sections for tt^ and K ± as functions of energy |30| to calculate the total number 
of nuclear interaction lengths traversed by a particle. The particle type, its energy and 
corrections to its momentum for energy loss through ionization were included. Given this 
information and the particle trajectory, we obtained the probability of punching through 
the absorbing material and producing track segments in the muon chamber. 

With the events in Fig. f|(a) that project to the CMU and CMP chambers, we assumed 
the third particle to be a pion and calculated its punch-through probability. We did similar 
calculations for K + and K~ . Using dE/dx information from the CTC, we determined that 
(56.0 ± 3.4)% of the third tracks are pions, where the uncertainty is purely statistical based 
on a fit. We assume charge symmetry for the relative numbers of K + and K~. The shapes of 
the mass histograms from all these calculations are nearly identical to each other and their 
sum is shown in Fig. §(a). The dominant contribution to the punch-through background is 
from K + because of its lower absorption cross section. 

As a check, we used this procedure to compute the number of K + and K~ punch-throughs 
from B — > J/ipK events and compared it with the actual number of punch-throughs in 
the data. For K + we predict 3.36 ± 0.46 events and observe 2 events. For K~ we predict 
0.65±0.08 events and observe 1 event. With such small samples, it is difficult to evaluate the 
systematic uncertainty and we arbitrarily assigned it a value comparable to these differences 
between the expected and observed number of J/ipK events. 

We estimate 0.88 ± 0.13 (stat.) ± 0.33 (syst.) events in the signal region due to hadron 
punch-through. 

2. Decay -in- Flight 

Pion or kaon decay- in-flight can contribute background to B c — > J/if> fi when a daughter 
muon from a meson decay is reconstructed as a track that projects to the J ftp decay point. 

We estimated this background from the events in the J/tp + track mass distribution 
shown in Fig. |(a). We assumed the third track to be a pion or a kaon and added it to a 
histogram with a weight that was the product of the following factors: 

• the probability that it would decay before reaching the muon chambers, 

• the probability that the data from the tracking system would be reconstructed as a 
track that points to the J/ip decay vertex. 

The decay probability is a simple calculation for each track. The probability for reconstruc- 
tion and vertex-pointing was calculated with a Monte Carlo program described in App. [A|. 
For the decay channels containing a J/i/j, the program forced pion or kaon daughters of a 
B to decay into a muon in the region upstream of the CMU chambers. It then traced the 
particles through the detector. This study included cases where the track did not originate 
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at the J/ip decay vertex, but decayed in a way that allowed a perturbed reconstruction 
which accidentally satisfied the vertex requirement. 

The events thus simulated were analyzed to determine the fraction of events for which 
the hadron and subsequent decay muon satisfied the muon identification criteria with a 
reconstructed track that projected to the J/ip decay point. The fraction depends only on 
the type of particle and on p?. The results of the calculation are shown in Figure [7] for 
kaons and pions. The ir/K ratio was determined from dE/dx as described in Sec. [IV A 1 . 
The appropriate fractions of the distributions for pions and kaons were added to yield the 
background mass distribution in Fig. |5|(b). 

The systematic uncertainty in the number of decay-in-flight background events arises 
from several sources: 

• uncertainties in the Monte Carlo calculation (12%), 

• uncertainties in the reconstruction efficiency for tracks from mesons that decay in the 
CTC (17%), 

• uncertainty in the tt/K ratio (10%). 

We estimate 5.5 ±0.5 (stat.) ±1.3 (syst.) events in the signal region due to the decay-in-flight 
background. 

3. Total False Muon Background 

The mass distributions for punch-through and decay-in-flight backgrounds are statis- 
tically indistinguishable in shape, and we have combined them for the likelihood analysis 
discussed in Sec. |V[ In the fitting region (3.35-11.0 GeV/c 2 ) we estimate 11.4 ± 2.4 (stat. 
© syst.) false muon events of which 6.4 ± 1.4 are in the signal region (4.0-6.0 GeV/c 2 ). 

B. False Electron Background 

Because of the requirement that the third lepton originates from the J/ip decay point, 
the main source of false-electron events among our B c candidates is B — > J/ijj + hadrons 
where one of the hadrons is misidentified as an electron. 

To determine the probability that a hadron was misidentified as an electron, we studied 
two independent sets of events deliberately chosen because they contain few real electrons: 
a dataset based on an inclusive jet trigger with a threshold transverse energy of 20 GeV 
(JET20) and minimum bias dataset based on a trigger that sampled beam crossings with 
no physics requirements (MB). 

The probability of misidentification of a track as an electron can depend on its transverse 
momentum and on the presence of nearby tracks. Therefore, we express this probability 
as a function of px and an isolation parameter J, defined to be the scalar sum of the 
momenta of particles within a cone AR < 0.2, divided by the momentum of the track under 
consideration. AR = a/ (Ar/) 2 + (A(f)) 2 is the radius of a cone in rj—cj) centered on that track. 
In this definition of isolation, a smaller / means more isolated. 
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The data in the JET20 and MB triggers contain a number of real electrons. In order 
to calculate the false electron probability for hadrons, the electrons were removed statisti- 
cally from the sample using dE/dx measurements. We computed the fraction f m of hadrons 
wrongly identified as electrons from the ratio of N e , the number of tracks satisfying all 
electron criteria, to N l , the number of tracks satisfying the purely geometric criteria. How- 
ever, a fraction f e of the tracks passing all electron criteria were, in fact, real electrons from 
heavy-flavor decays and from conversions, i.e., pair production by photons and Dalitz pairs 
as discussed in Sec. [IV C| . From dE/dx measurements we found f e to be 0.74 ± 0.02 in the 
JET20 data and 0.64 ± 0.07 in the MB data. Thus, 

/m=^x(l-/ e ). (5) 

Figure shows f m as a function of px for the two data sets and for two ranges of the isolation 
parameter. The results from the MB data differ from those of the JET20 data by 10%, and 
we adopted this as a measure of the systematic uncertainty in this calculation. 

We calculated the number of background events due to misidentified hadrons in J/ipe 
(Fig. 0(b)) by selecting J/ip + track events (Fig. 0(a)) in which the third track is required 
to satisfy the purely geometric criteria for electron identification. For each such track, 
we calculated / and weighted its contribution by the probability f m (pT,I) determined in 
the JET20 studies. A mass histogram of the weighted sum is given in Fig. 0(a). The 
number of hadronic background events determined with this technique was consistent with 
that expected from the dE/dx distribution data prior to the application of the dE/dx 
requirement. Figure |10|(a) shows the results of a dE/dx calculation applied to the third 
track for events in Fig. 0(a). Most are hadrons. These tracks were then required to satisfy 
all the electron identification criteria except the dE/dx requirement. Results of the dE/dx 
calculation for the surviving events are shown in Fig. |i~0|(b). For most of the surviving events, 
the third track is an electron. 

We estimate 2.6 ± 0.3 (stat. © syst.) events in the signal region due to false electrons 
and 4.2 ± 0.4 such events in the fitting region. 

C. Conversion Background 

Photon pair production in material around the beam and Dalitz decays both produce 
e + e~ pairs. The reconstructed track for one member of a pair can pass through the J/i/j 
decay point and be selected as a candidate for B — > J/ipX. After applying other electron 
identification criteria and the vertex constraint (Sec. plf) , we found and rejected two such 
"conversion" events by searching for the partner track in the J/ip + track sample with 
ct* > 60 /im. However, a track can contribute to the background in the J/ipe events if its 
partner track has low momentum and escapes detection. 

To estimate the magnitude and shape of this background in the M(J/ipe) distribution, 
we performed a hybrid Monte Carlo calculation based on the J/ijj + track events. The Monte 
Carlo program replaced the third track in the event by a 7T°. It forced 1.2% of the 7r°'s to 
decay through the Dalitz channel and the rest through two-photon final states. The program 
propagated the photons through the surrounding material with tabulated probabilities for 
e + e" production, and it propagated the resultant charged particles through the detector 
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simulation. We used each event 100 times, rotating its azimuth by a random angle to 
sample all parts of the detector. Figure 11 shows the momentum spectrum for the track 
which fulfilled the requirements for the third lepton and the spectrum for the other member 
of the pair. These hybrid events were subjected to the B c analysis procedures. Roughly half, 
(48.6 ± 1.9)%, of these "conversion" background events were rejected because the partner 
was detected. Thus the ratio of undetected or residual conversions to detected conversions 
is R ce = 1.06 ±0.08 (stat.). 

In the simulation, the J/ipe mass distributions arising from detected and undetected 
conversions have the same shape. Fig. [|(b) shows this shape normalized to an area equal to 
the expected 2.1 undetected conversion background events. 

Systematic uncertainties arise from statistical uncertainty in the efficiency for finding 
the conversion partner (28%), from uncertainty in the shape of the J /if} + track mass 
distribution for these events (9%), and from differences in p T distributions between the data 
and the sample used to calculate this background (13%). Combined, they are 32%. 

The statistical uncertainty from two events is the largest contribution to the overall 
uncertainty in the conversion background, and we quote the Gaussian approximation of the 
uncertainty here. In the likelihood analysis of Sec. |V|, the two detected events, N ce = 2, enter 
as a Poisson term. The systematic uncertainties are incorporated in the ratio of undetected 
to detected conversions, R ce = 1.06 ± 0.36 (stat. © syst.). The residual background is the 
product N ce R cv . 

The mass distribution for the conversion background distribution in Fig. |^(b) contains 
2.1 ± 1.7 (stat. © syst.) events in the fitting region. Of these 1.2 ± 0.9 events are in the 
signal region. 



D. BB Background 

BB pairs produced during the pp collision can mimic the B c — > J/vpiv signature when 
a B decays into a J/if> and its associated B or any of its daughters decays into a lepton. 
If the lepton track projects through the J/ip vertex, the event may not be distinguishable 
from a B c decay and would be a part of the irreducible background. 

The BB background was determined by a Monte Carlo simulation (App. |A]). One B 
was required to decay into a final state containing a J/ip, and the other B was allowed 
to decay through all channels. We simulated the detector response, and we required the 
simulated events to pass the di-muon trigger criteria. To avoid double-counting false-lepton 
backgrounds, we eliminated candidates where the third track was a hadron. We then per- 
formed the B c analysis on these events. We used B — > J/if> K events to normalize the Monte 
Carlo simulation to the data. The resulting mass distributions are shown in Figs. ||(c) and 
1(c). 

The systematic uncertainties in the estimate of this background include the trigger simu- 
lation (5%), the uncertainty in the branching ratio B + — > J/if> K + (10%), and Monte Carlo 
statistics (11%). 

We estimate 0.7 ± 0.3 (stat. © syst.) J/ip /i events and 1.2 ± 0.5 J/ipe events in the 
signal region due to BB background. The corresponding numbers in the fitting region are 
1.44 ± 0.25 J ftp fj, events and 2.3 ± 0.9 J/ip e events. 
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E. Other Backgrounds 



We have considered three additional potential sources of background to the decay B c — ► 
J /if) I v. They are 

• false J /if) candidates from the continuum background of the di-lepton spectrum, 

• J /if) + cc production in which the charm decays semileptonically, and 

• decays of as yet undiscovered baryonic bc states such as the E bc . 

We estimate that these make negligible contribution to our background. 

The false J / if) background is very small after mass and vertex constraints are applied to 
the data. We selected two side bands in the J /if) mass distribution. In each we substituted 
the central mass for the side band in our fitting procedures. We found 3 u J/ifj" + track 
events that satisfied our track selection criteria. In none of these did the third track satisfy 
our criteria for muons or electrons. The dominant source of false J / if) candidates is B decay 
to a real muon along with a hadron falsely identified as a muon because of punch-through 
or decay- in-flight (Sec. |LV A ). Either the associated B or a daughter D has a branching 



fraction of roughly 0.1 for yielding a third lepton. The probability for another hadron 
falsely identified as the third lepton is even lower, roughly 0.01. Our background estimate 
is 3 x 0.1 x 0.5 = 0.15, where the factor 0.5 is the ratio of widths for the central peak vs. 
two side bands. The 90% confidence upper limit on 3 events is 6.7 events, which yields an 
upper limit of 0.34 events. We neglect this source of background. 

It is possible for additional charm to be produced along with prompt J /if) mesons with 
production mechanisms similar to those for B c production. Several of our selection require- 
ments suppress background from such events in which the additional charm decays semilep- 
tonically. As is the case with the BB background, the prompt J /if) + charm background 
is suppressed because the J /if) and lepton do not generally form a common vertex. Addi- 
tional suppression of charm-daughter leptons results from the isolation cut in the electron 
channel and the high transverse momentum requirements in both channels. Finally, since 
these events are prompt, they mostly fail the ct* requirement. For the lifetime meaurement 
discussed in Sec. [VI I| , we studied the ct* dependence of the signal and various backgrounds. 



They account for the distribution of candidate events at low ct*, and there is no evidence 
for additional background from J /if) + charm. Therefore, we neglect it. 

The as yet undiscovered hyperon E bc can decay into a tri-lepton topology, e.g. E bc — > 
J /ip S c followed by S c — > H I v. The production cross section for such a particle is likely 
to be significantly less than that for the B c . Alternate standard-model decay modes for 
S& c fail our B c identification criteria. The same observation can be made for other baryons 
containing a b quark. We assumed no background from these particles. 
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V. THE MAGNITUDE AND SIGNIFICANCE OF THE B c SIGNAL 



Monte Carlo calculations (App. for B c — > J/ipfi and B c — ► J /ip e with the B c mass 
assumed to be 6.27 GeV/c 2 yielded the tri-lepton mass distribution shown in Fig. |l2"|(a). 
The normalization anticipates the results of the fit described below. The electron and muon 
mass distributions are used seperately, but the figure shows the combined distribution since 
the differences are small. We assume equal branching fractions for the two decay modes, and 
we expect the ratio of J/ipe to total B c — > J/ip I v events to be given by the efficiency ratio 



R £ = 0.58 ±0.04 discussed in Sec. [11 fD| The mass distribution for the sum of the normalized 
backgrounds for muons and electrons is shown in Fig. 0(b). The mass distribution for all 
B c candidates is shown in Fig. 0(c). 

The expected background is unable to account for the observed data distribution. In 
order to test this statistically and to determine the magnitude of the signal needed to account 
for the excess, we adopted two approaches. The first was a simple "counting experiment" 
based on the number of events in the J/ip + lepton mass range from 4.0 to 6.0 GeV/c 2 . 
However, this ignores additional information in the shapes of the distributions and the 
yield in the extended mass range populated by backgrounds but not by signal. Our second 
approach employed a binned likelihood fitting procedure that includes the shape of the 
distributions over the full mass range, 3.35 to 11.0 GeV/c 2 . To account for the excess in 
the data over expected background, the fit varied the normalization of the signal shape of 
Fig. |T2"Ka) and calculated its uncertainty. The bins are those shown in Figs. |3| and || except 
that the lowest bin in the figures, 3.05 to 3.35 GeV/c 2 , was not used in the fit. 

In both approaches, we computed the probability that a random fluctuation of the back- 
ground is sufficient to account for the observed data in the absence of a B c contribution. 
This is the "null hypothesis." 

We also performed an unbinned likelihood analysis using spline fits to the parent distri- 
butions. The results are completely consistent with the binned likelihood analysis. We also 
varied the assumed B c mass from 5.28 to 7.52 GeV/c 2 . Within the range 6.1-6.5 GeV/c 2 , 
which embraces all the theoretical predictions, we found the fitted number of B c events to 
be insensitive to the assumed mass These issues are discussed in Sec. |VT[ 



A. The Counting Experiment 

In the signal region of J/ipt mass, we observe 19 J/ipe candidates and 12 J/ip /i candi- 
dates. Table | summarizes the backgrounds from the various sources of background discussed 
in Sec. [TV]. The expected total backgrounds are 5.0 ±1.1 events for J/ip e and 7.1 ± 1.5 events 
for J/ip /i, leading to a combined signal of 18.9 ± 5.6 events. From these results, we tested 
the null hypothesis by folding the Gaussian uncertainties in the estimated mean number of 
background counts with their Poisson fluctuations. This allowed us to determine the prob- 
ability that the backgound would fluctuate up to the observed number of events. The null 
hypothesis probabilities are 2.1 x 10~ 5 for the J/ipe sample and 0.084 for the J/ip fi sample. 
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B. Likelihood Analysis: Fit to the B c Signal 

We used a normalized log-likelihood function for testing and fitting our data and back- 
ground estimates. It used the shapes of the distributions over the mass range 3.35 to 11.0 
GeV/c 2 , and it included as input all the information on the tri-lepton mass distributions 
for signal and for background discussed in earlier sections. The likelihood function has a 
necessary and sufficient set of parameters to fit these distributions to the observed data. It 
also included constraints such as the expected fractions of events in the two decay channels. 

In Appendix |C| we discuss the normalized log-likelihood function £ 2 = — 21og(£/£o) 
used to fit our data, where £ is the likelihood function and Co is its value for a perfect fit. 
Maximum likelihood is equivalent to minimum £ 2 which has properties similar to those of 
X 2 '■ The only unconstrained parameter in the fit is n' e , the total number B c — > J/tpn and 
B c —* J/ipe events in the fitting region, i.e. in the J/ip£ mass range 3.35-11.0 GeV/c 2 . All 
other parameters in the fit are constrained by externally derived information. 

At the minimum in £ 2 , the number of J/ipl events in the fitting region is 

n" = (n'" + n' e ) = 20.4+| 2 events (6) 

with E^/Nd.o.f. = 38.1/26, where Nd.o.f, is the number of degrees of freedom in the fit. In 
the Monte Carlo signal distribution in Fig. [1^, (93.0 ± 0.6)% of the events fall in the signal 
region (4.0-6.0 GeV/c 2 ). We scale 20.4 events by this value to calculate 

n e = (n^ + n e ) = 19.0+^ events (7) 

in the signal region. This is in excellent agreement with the counting experiment result. 

Figure shows the contributions to the background and signal for B c — > J/ipeX and 
B c —>■ J/ip fiX separately resulting from the binned likelihood fit, and Fig. |14] shows the 
combined data. 

Figure [15| shows £ 2 plotted as a function of the assumed number of B c mesons in the 
data sample. For each value of n' e , £ 2 was minimized as a function of the other parameters. 
Table shows the input constraints and fitted values for the background normalizations 
and for other parameters. 

To evaluate the quality of the fit, we observe that, to the extent that £ 2 behaves like \ 2 i 
P(£ 2 ) = 8%. We made a more reliable estimate of this probability by generating a large 
number of Monte Carlo "pseudo-experiments." First, we generated random backgrounds 
with Gaussian-distributed uncertainties based on the shapes and normalizations determined 
in Sec. [TV[ To this we added a signal contribution with the fitted magnitude varied according 
to the uncertainty from the fit. Bin- by-bin, the signal plus background value served as the 
mean for a number of events randomly generated according to a Poisson distribution. This 
constituted a pseudo-experiment with a B c signal. We ran the fitting program on each 
pseudo-experiment. The £ 2 distribution for these is shown in Fig. |16[ The probability of 
finding £ 2 > 38.1 is 5.9%. 

Only two assumptions about the B c signal distribution were used in the fit: the B c 
mass and the relative contributions to the electron and muon channels. The choice of 6.27 
GeV/c 2 for the mass will be considered in Sec. VI. As a test we fit the data with the electron 



fraction r £ allowed to vary freely, not constrained to R e = 0.583 ± 0.043. The results of this 
fit were: £ 2 = 37.7; the number of signal events was 20.3, and the fitted electron fraction 
was r £ = 0.65 ± 0.14, consistent with R £ . 
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C. Likelihood Analysis: The Null Hypothesis 



The null hypothesis is the postulate that there is no B c signal and that a statistical 
fluctuation in the background is responsible for the apparent excess in the data. In order 
to test this, we again computed the results for a large number of "pseudo-experiments" or 
trials in the manner described above, except that we omitted the signal contribution. With 
n' e allowed to vary, we ran the fitting program to return the fitted number of B c events in 
a distribution devoid of real signal. Figure [T7| shows a histogram of n' e for 351,950 pseudo- 
experiments. The fitted signal tends to compensate for statistical fluctuations, positive or 
negative, from the correct background shape. The peak at zero events includes those trials 
consistent with a negative contribution from the signal distribution. No pseudo-experiments 
gave values of n' e exceeding 20.4. We extrapolated the fitted shape of the distribution and 
estimate its area above 20.4 to be 0.22+°;^° out of 351,950 trials. Thus, the probability that 
a random fluctuation of the background could produce the observed data distribution is 
0.22/351,950 = 0.63 x 10~ 6 . This is equivalent to 4.8 standard deviations in significance. 

In the following sections, we assume that the excess events are due to the existence of the 
B c meson. We describe measurements of its mass, its lifetime and its relative cross section 
times branching fraction, all of which are consistent with values expected for the B c . 
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VI. THE B c MASS 



In order to check the stability of the B c signal, we varied the value assumed for the B c 
mass. With the procedures described in Sec. |V| and App. |A], we generated Monte Carlo 
samples of B c — > J/ip i v with various values of M(B C ) from 5.52 to 7.52 GeV/c 2 . For each 
of these samples, we propagated the B c — > J/ijjlv final-state particles through the detector 
simulation programs to obtain the tri-lepton mass spectrum, i.e. a signal template. The 
signal template for each value of M(B C ) together with the background mass distributions 
was used to fit the mass spectrum for the data. The best-fit log-likelihood value shows a 
rough parabolic dependence on the assumed B c mass, and this yields a measurement of 
M(B C ). 

We performed this analysis with the binned log-likelihood analysis described in Sec. [V] 
and with an un-binned log-likelihood analysis. The two methods yielded nearly identical 
results, but the binned method exhibited slightly more scatter about a smooth dependence 
on mass. We present the unbinned results here because this method is not sensitive to 
binning fluctuations. 

For each assumed B c mass, a signal template was formed with a smooth spline fit to the 
Monte Carlo distribution. Figure [18| shows the generated distributions and spline fits for a 
sample of the templates used in this study. Background templates formed in the same way 
were independent of the assumed B c mass. Most contributions to the unbinned log-likelihood 
function were the same as those in Sec. [V] and App. |C 2| for the binned fit. However, the sum 
over bins of Poisson terms was replaced by of sum over events of log-probabilities. This is 
discussed in App. |C 3| . In this analysis we compare the log-likelihood to its value at minimum 
C m in, and we define the relative log-likelihood function £^ as a function of M(B C ). 

& = -2]n(JL-). (8) 



C 



mtn 



At each assumed value of M(B C ), several Monte Carlo samples and corresponding signal 
templates were generated in order to determine the sensitivity of the fit to statistical fluc- 
tuations in the Monte Carlo simulation. This provides us with an uncertainty on the values 

Of Cr 

Figure 0(a) shows the dependence of £^ on M(B C ). The figure includes a parabolic fit 
to The parabolic fit yields a best fit value of 6.40 GeV/c 2 with a statistical uncertainty 
of 0.39 GeV/c 2 . 

As in Sec. |V| we generated a sample of pseudo-experiments based on the fitted results 
with the assumed mass of 6.27 GeV/c 2 . The distributions of M(B C ), its uncertainty, the 
number of B c and its uncertainty were consistent with the results in the experimental data. 
This provides some confidence that the model used to fit the data is adequate to the task. 
The comparison between the unbinned log-likelihood function for the experimental data 
and that for the pseudo-experiments was closely similar in shape and width to that for the 



binned likelihood analysis (Fig. [16]) 



We considered a number of sources of systematic uncertainty in this measurement: 

• distortion of the signal mass distribution arising from decay to higher-mass cc states 
rather than J/ip (0.09 GeV/c 2 ). 
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• fitting procedures, estimated from the difference between binned and unbinned anal- 
yses (0.08 GeV/c 2 ), 

• finite Monte Carlo statistics in the signal template (0.04 GeV/c 2 ), 

• variations in the B c mass distribution due to 6-quark production spectrum (0.02 
GeV/c 2 ), 

• Monte Carlo simulation of the CDF trigger (0.02 GeV/c 2 ), 

These uncertainties are small in comparison with the statistical uncertainty. In quadrature, 
they sum to 0.13 GeV/c 2 . 

Figure |l9|(b) shows that the magnitude of the B c signal is stable over the range 
of theoretical predictions for M(B C ), and our experimental measurement of the mass is 
M(B C ) = 6.40 ± 0.39 (stat.) ± 0.13 (syst.) GeV/c 2 . 
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VII. THE B c LIFETIME 



We extended our analysis to obtain a best estimate of the mean proper decay length cr 
and hence the lifetime r of the B c meson. The information to do this is contained in the 
distribution of ct* which is defined in Eq. 0. We changed the threshold requirement on ct* 
from ct* > 60 fim to ct* > —100 /im. This yielded a sample of 71 events, 42 J/ipe and 29 
J/i/j [i. We determined a functional form for the shapes in ct* for each of the backgrounds 
(Fig. ^1|). To these, we added a resolution-smeared exponential decay distribution for a B c 
contribution, parametrized by its mean decay length cr. Finally, we incorporated the data 
from each of the candidate events in an unbinned likelihood fit to determine the best-fit 
value of cr. 

Since the neutrino in B c —>■ J ftp iv carries away undetected momentum, the true proper 
time for the decay of each event cannot be calculated from ct* . The relationship between 
ct* and ct is: 



ct' = I (9) 



where K for an event is given by 



M(B C ) pr{JH>t) m 

We assume M(B C ) = 6.27 GeV/c 2 , but pt{B c ) is unknown for single events, and therefore, 
we cannot correct for K event-by-event. In an ideal data sample with no background and 
a known pt(B c ) distribution, one finds (ct*) = (ct)(l/K) = cr(l/K), where (ct*) is the 
average over the data, and (l/K) is the average over pt(B c ) and pr(J W -0- 

For B c — > J/ipe and B c — > J/ipfi, we obtained the K distributions H(K) by Monte 
Carlo methods. Figure ^0] shows the results of these calculations for the kinematic criteria 
pr(e) > 2 GeV/c or p T (/i) > 3 GeV/c, and 4 GeV/c 2 < M(J/ip£) < 6 GeV/c 2 . Since the 
criteria differ for the electron and muon, the fT-factor distributions for these channels were 
determined separately. For the exponential dependence of Rk on (1/cr) (Sec. pi D|) , the 
distributions in Fig. |20] can be adequately represented by (l/K) = 0.88 ± 0.02, where we 
have adopted the difference between the two distributions as the uncertainty. 

The quantity ct* was determined for each event by the relation given in Eq. [I]. The points 
with uncertainties in Fig. E2I show the binned ct* distributions for the J/ipe and J/ipfi data. 



The two decay channels are combined in Fig. 23 



A. Background and Signal Distributions in ct* 

We used a procedure similar to that described in detail in Ref. |^9[ to account for back- 
grounds. We constructed functions to represent the ct* distributions, for signal and back- 
grounds and convoluted them with a Gaussian resolution function. 

The evaluation of backgrounds for events with ct* greater than 60 /im was described in 
Sec. [TVl The same procedures were used independently for events with ct* between -100 fim 
and 60 /im which have "prompt" contributions from direct charmonium production. 
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We obtained the best fit to the ct* distributions for each of the backgrounds using the 



same methods discussed in Sec. [IV| for the background rate determinations. The general 



shape in x = ct* used for each of the backgrounds was a sum of three terms: 

• A right-side (ct*> 0) exponential dominated by the decay of ordinary Bs in the back- 
ground. Its fractional contribution is /+ and its exponential slope is A+. 

• A left-side (ct*< 0) exponential to account for an observed low level background from 
daughters of B decay incorrectly associated with particles from the primary intraction 
vertex. Its fractional contribution is f J _ and its exponential slope is X 3 _. 

• A central Gaussian to account for prompt decays. Its fractional contribution is (1 — 

r. •-/'!. 

The index j stands for the various background contributions from false muons (j = ffi), false 
electrons (j = fe) and undetected conversion electrons (j = ce). For the BB backgrounds 
(j = Be), the central Gaussian term in Eq. |ll| was not needed, i.e. f+ e + f^ e = 1. The 
exponentials were convoluted with a Gaussian resolution function. This sum can be written 



F(x) = (1- f> + - fL)G(x;s>a 
A 




+^j-9(x) exp 
+ 

+j=-6(-x) exp ( ) <g> G(x; sV). (11) 

where the Heaviside function 9{x) is defined as 8(x) = 1 for x > and 8(x) = for x 

< 0. The product s^a is the one-standard-deviation width of the Gaussian distribution, 

where a is the measurement uncertainty on x for each event and is a fitted scale factor. 
In all background fits, the s J were consistent with a common value of s = 1.4. Therefore, 



s was fixed at that value. Figure |21] shows the distributions and fitted functions for the 
backgrounds. Table [III] shows the fitted shape parameters for each background. The values 
of A^ suggest that the backgrounds are dominated by partially reconstructed B mesons. 
Table |I| also shows the numbers of events for each background. These differ from the 
corresponding numbers in Tables | and [TI| because of differences in the selection criteria for 
ct* and tri-lepton mass used here. For this reason, we adopt a double-prime notation for 
this analysis, e.g. n"^ for the number of false muon events with M(J/ip£) in the range 4.0 
to 6.0 GeV/c 2 and with -100 //m < ct* < 1500 [im. 

Our fitting procedure accounted for a difference between the relative pion and kaon 
fractions contributing to the prompt background and that contributing to background in 
the 5-like region with ct* > 60 /im. The fit also allowed variation in the relative probability 
for pions and kaons to be falsely identified as electrons or muons. These considerations allow 
additional variation of the values of f± in Table |II| and are discussed in App. |C 4| . 



We assumed an exponential decay for the contribution from B c , but we convoluted it with 
the K distribution and a Gaussian distribution to account for measurement uncertainty. 



F e S ig{ x ' CT ) 



H(K) (Jpj e("^) ® G(x; s l a) 



dK (12) 
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where £ = //, e. The weighted sums of signal and background probability distributions are 
defined in App. |C 4| . 



B. Unbinned Likelihood Fit for ct 

We used an unbinned likelihood method to obtain a best estimate of cr for each decay 
channel individually and for the combined dataset. A parameter in the fit was assigned 



to each of the quantities in Table |T|. The numbers of events in each background were 
constrained by their measured or calculated values as in the previous sections. The full 
covariance matrices from the fits that determined the background shape parameters were 
used to constrain them in the lifetime fit. As before, we used the total number of events n" 1 
and the electron fraction r e to describe the B c signal with n" e = r e n" e and n"^ = (1 — r e )n' . 
The only parameter unconstrained by information beyond the candidate events was cr, the 
mean decay length for the B c contribution to the ct* distribution. The likelihood function 
is presented in App. 



The result of the log-likelihood fit to the ct* distribution for J/if)e events is 

CT = 122±%tim (13) 

For J '/if) \i events, the fit yielded 

cr = 1 72+*° fim (14) 
The solution for a simultaneous fit to all events is 

CT = 137+53 (15) 



0.46±g;il ps (16) 



The variation of —2 ln(£) from its minimum as a function of cr is shown in Fig. |24|. The 
simultaneous fit also determined the number of B c events to be 

n" 1 = 34. 2±f;§ events (17) 

With the mean decay length above, the acceptance for ct* greater than 60 fxm is 0.611q; 15 , 
and we can calculate 

n l = 20.9±H events (18) 

for comparison with Eq. [7[ Clearly there is a large correlation between these two numbers 
because of the largely overlapping event samples. However, the consistency of the size of the 
B c signal as determined from both the tri-lepton mass distribution and the ct* distribution 
adds confidence to the result. 
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C. Statistical Tests of the Fit 



In order to test the adequacy of our model for signal and background, we ran a number 
of pseudo-experiments based on the fitted values of R e , n , and the background parameters. 
For each of the pseudo-experiments, we varied these parameters randomly according to the 
appropriate Poisson or Gaussian uncertainties. The value of cr was fixed at 140 /im for all 
pseudo-experiments. From these quantities, we constructed the J/ipe and J '/if) /i probability 
distributions for the independent variable ct*. The dataset for a pseudo-experiment consisted 
of contributions from a signal plus three types of background for J/if>e and a signal plus two 
types of background for J /if) /i. For each of the five backgrounds the number of events was 
allowed to fluctuate according to Poisson statistics, and the value ct* was chosen randomly 
according to the appropriate probability distribution. The total number of signal events 
was chosen according to Poisson statistics, and each event was designated J /if) e or J /if) fj, 
with probability determined by R e . These samples were then subjected to the same fitting 
procedures as the experimental data. The comparison between the results for the pseudo- 
experiments and those for the data tests the adequacy of the fitting function to represent 
the data. 

Figure ^5|(a) shows the distribution for the log-likelihood with a mean value of —382 
and an r.m.s. width of 49. The experiment yielded —430, which corresponds to an 84% 
confidence level. Figure p5|(b) shows the distribution of fitted values of cr. The mean of the 
distribution, 144 /im, agrees closely with the input value of 140 /im, and the width is 44 /xm, 
which consistent with the measured uncertainty. Figure 0(c) shows the distributions of the 
upper (solid histogram) and lower (dashed histogram) uncertainties from the fits. Arrows 
indicate the corresponding uncertainties from the experimental data. They are in reasonable 
agreement with the results from the pseudo-experiments. Figure |25](d) shows the distribution 
for deviation of the fitted cr from the input value normalized to the uncertainty from each 
fit. 

We conclude that the model used to fit the data is adequate and that the resulting 
log-likelihood value and fitting uncertainties are consistent with expectations based on the 
uncertainties in the data. 

D. Systematic Uncertainties 

The uncertainty reported by our fitting program already includes some sources of sys- 
tematic uncertainty because of the way we constrained the parameters describing the signal 
and backgrounds. The fit shows a correlation of —10% between cr and the prompt elec- 
tron fraction discussed in App. |U4| . The correlations with all other fitting parameters are 
less than 5%. Thus, the cr value varies only a fraction of a standard deviation as other 
parameters in the analysis are varied. Refitting with parameters fixed at values different 
from nominal gives results consistent with this. We estimate the systematic uncertainty 
included in the fitting uncertainty to be less than 10 /im. Thus, the fitting uncertainty is 
overwhelmingly statistical, and we quote it as such. 

Below we discuss additional sources of systematic uncertainty. Combined in quadrature, 
they amount to about one-fifth the statistical uncertainty. 
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The K distribution (Eq. |10| and Fig. |20|), which was used to compensate for the infor- 
mation lost by our inability to detect the neutrino, is vulnerable to errors in our model of 
the B c production spectrum and its decay kinematics. 

Figure [| shows that the pt spectra for data and background are very similar to that 
calculated for B c which was used to generate the K distribution. To generate the B c Monte 



Carlo events, we used the next-to-leading order calculation of the b quark spectrum |3l|j32 



with the MRSDO parton distribution functions (PDF) |33|], m;, = 4.75 GeV, and the renor- 



malization scale ji = /xq = \/ m & + Pr(b) 2 . We also generated a B c Monte Carlo sample using 



the CTEQ4M PDFs [34] to obtain a new K distribution and used it to fit the signal sample. 
The value of cr thus obtained differed by 2 /im from the value in Eq. [15|. Therefore, we 
assign ±2 /im systematic uncertainty for the PDFs. 

We also refit the data with the assumed B c mass changed by ±150 MeV. This yielded a 
variation in cr of ±1.6 /im. 

A B c can decay to a lepton, a neutrino, and a higher mass cc state that can subsequently 
decay to J/ipX. This would satisfy the requirements for a candidate event, but would give 
rise to a different K distribution. Calculations based on the ISGW model |3J indicate that 
the largest such contribution comes from B c — » ip(2S)£v, which could account for 12% of the 
B c — ► J/ip £ v candidate sample. We generated events of this type to obtain a K distribution 
that we used to refit the candidate events. The value of cr changed by 1.9 /im which we 
adopt as a measure of the systematic uncertainty for this effect. We also considered the 
effects of B c — > J/iprv , B c — > J/ip D s , and B c — > J/ipD*. We estimate their contribution 
to the B c sample to be less than 5%. We assume that they produce no change in the lifetime. 

Our model for B c decay |35j uses a V — A matrix element. As alternative, we generated 



events with the ISGW model |36| to obtain a new K distribution and refit the data. This 
indicates a possible systematic uncertainty of ±2.0 /im 

In order to test possible bias in our experimental trigger, we turned off the trigger 
simulation in our Monte Carlo program and generated a sample of events without it to 
obtain a K distribution. We assign ±1 /im uncertainty for this effect. 

For each event in the lifetime analysis, the raw uncertainty in ct* was multiplied by a 
scale factor, s — 1.4 that best fits the distributions in our background studies. We changed 
this factor by ±0.4 and re-fit the background shapes. We assign a systematic uncertainty 
in cr of ±8.3 fim for this effect. 



In another analysis of B hadron lifetimes |29 ], we studied the effects of detector alignment. 



From this work, we assign an uncertainty on cr of ±2.0 /im. 

In quadrature, these uncertainties sum to ±9.4 /im, and we quote this as our systematic 
uncertainty with the caveat that some other sources have already been included in the fitting 
uncertainty which, nevertheless, remains predominantly statistical. Thus, our result is: 

cr = 137^49 ( stat -) ± 9 (syst) V m ( 19 ) 
r = 0.46+g;i| (stat.) ± 0.03 (syst.) ps (20) 
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VIII. B c PRODUCTION 



From the event yield of Sec. [V], we calculated the B c production cross section times 
the B c — > J/if)lv branching fraction a ■ BR(B+ — > J/ip£ + u). We express this product 
relative to that for the topologically similar decay B — > J /if) K because the systematic 
uncertainties arising from the luminosity, from the J /if) trigger efficiency, and from the CTC 
track-finding efficiency cancel in the ratio. Our Monte Carlo calculations yielded the values 
for the efficiencies that do not cancel in the ratio. We assumed that the branching fraction 
is the same for B c — > J /if) e and B c — > J /if) /i. 

We use the number of B c events from Eq. |6] and the number of J/if)K events from the 
fit in Fig. |. 

n» + n e = 20.4^-2 events (21) 
n K = 290 ± 19 events (22) 



In order to be consistent with the efficiency calculations of Sec. |III D| , the B c event count 
is that for M(J/ip£) in the range 3.35 to 11.0 GeV/c 2 . We relate these quantities to the 
luminosity £, to the products of cross section and branching fraction a ■ BR, and to the 
efficiencies discussed in Sec. [Ill D| . 

n e = C ■ a{B c ) ■ BR(B C -> J/ip tv) ■ e e (23) 
= £ ■ a(B c ) ■ BR(B C -> J/ip tv) ■ (24) 
K - C ■ a{B) ■ BR(B -> J/if) K) ■ e K (25) 



n 



n e + n^ _ a(B c ) ■ BR(B C -> J/ip tv) e e + 

n K ~ a(B) ■ BR(B -»• J/if)K) ^ ' 

a{B c ) ■ BR(B C -> J / if) tv) R K 



a(B) ■ BR(B — >• J/if)K) R e 



(27) 



We used the value of R e from Eq. [| We calculated the efficiency ratio R K from Eq. | and 
the lifetime discussed in Sec. |VII| to be 



R K = 0.263 ± 0.035 (syst.) (lifetime). (28) 



As was discussed in Sec. VII, there can a contribution to our data sample from other decay 



modes of the B c . Estimates of partial widths for higher charmonium states |!6| yield an 
upper limit of 12% for their contribution to the signal. The estimated contributions from 
final states involving D s , D*, and r with subsequent decay to e or fi total less than 5%. We 
assume an uncertainty equal to the magnitude of the correction 1/1.17 = 0.85 ± 0.15. With 
these values we find 

«(^ t ).f'' M ^f': ) (29a) 
y 1 a{B) ■ BR(B — > J/ipK) v ; 

n e + n M R £ . , „ 

= 0.85 — ■ — 29b 

n K R K y ' 

= 0.132+°;°^ (stat.) ± 0.031 (syst.) toml (lifetime). (29c) 



30 



The statistical uncertainty is from the event counts and the systematic uncertainty is from 
the efficiency ratios and the correction for other decay modes. 

Based on Monte Carlo studies, the effective kinematic limits for B c mesons in this study 
are: transverse momenta pt > 6.0 GeV/c and rapidity \y\ < 1.0. 

Figure |26] shows theoretical predictions of the ratio TZ{J/ipiv) as a function of the as- 
sumed lifetime of the B c . The shaded regions in the figure represents the prediction and 
its uncertainty for two different assumptions about the semi-leptonic width T s j_ = T(B+ — > 
J/ipl + v). Assumed in the theoretical predictions are 



V{Bt 
orT(.B+ 



v bc 

J/ip£ + v) 
a(B+) 



BR(B~ 



a(b) 
a(B + 

J/i)K + 



0.041 ±0.005 [30], 


(30) 


(30.6 ± 16) x 10~ 15 GeV 0, 


(31) 


16.5 x 10~ 15 GeV @, 


(32) 


1.3 x 10~ 3 13j, 


(33) 


0.378 ±0.022 [30|, 


(34) 


(1.01 ±0.14) x 10 -3 0. 


(35) 



Fig. ^6] also shows the measured cross section ratio (Eq. |29c|) plotted at the measured value 
of the lifetime. 

In Sec. | we referred to results from previous searches for the B c meson through its decay 
to various final states (f.s.) including J/ipir, J/ip7i + 7i^7i, J/ipai and J/ipiu. We have 
converted the upper limits quoted in these searches to calculate in each case a corresponding 
upper limit on TZ(f.s.) as defined in Eq. |29a| . For these conversions, we used BR(Z — > 
bb) = 0.1546 ± 0.0014, BR(Z -> qq) = 0.6990 ± 0.0015, BR(b -> B + ) = 0.378 ± 0.022, 

30| . The limits reported for the LEP experiments 



BR(B + -> J/ij) K~ 



;i.01±0.14)xl0" 



are for the sums of the two charged conjugate modes, and they are modified by a factor of 



2 for this calculation. Table IV shows the results of these calculations. 
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IX. SUMMARY AND CONCLUSIONS 



This paper reports the observation of B c mesons. The decay mode used for the study 
was B c — > J 'ftp IX where I is either an electron or a muon. A total of 31 events for which the 
mass of J/ip£ system was between 4.0 and 6.0 GeV/c 2 were found. We performed a detailed 
study of backgrounds and estimate their contribution to this sample to be 12.1 ±1.9 events. 
In the wider mass range 3.35 to 11.0 GeV/c 2 we found 37 B c candidates with an estimated 
background of 21.4 ±3.1 events. We performed a shape-dependent likelihood fit to the mass 
distribution and found that it required a B c contribution of 20Atfl of which 19.0jl|;f have 
masses between 4.0 and 6.0 GeV/c 2 . A fit without a B c contribution was rejected at the 
level of 4.8 standard deviations. 

By repeating the above procedure with a number of assumed masses between 5.52 
GeV/c 2 and 7.52 GeV/c 2 we determined that the mass of the B c meson is M(B C ) = 
6.40 ± 0.39 (stat.) ± 0.13 (syst.) GeV/c 2 . 

We studied the displacement of the B c decay vertex position from the average beam line, 
and from it we measured the B c lifetime to be t(B c ) = 0.46+Q'Jg (stat.) ± 0.03 (syst.) ps . 

Finally, we estimated ratio of the product of the production cross section times branching 
fraction for B+ — > J/ipi + v to that for B + — > J/ipK + to be 

a(B+).BR(B+^J/il>Pv) ft1oo+0 . M1 ,^ ^ n nQ1 +0.032 



a(B+) ■ BR(B+ -> J/ipK+) 



0.1321^ (stat.) ± 0.031 (syst.) (lifetime). 
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APPENDIX A: EVENT SIMULATION 



A number of quantities and distributions needed for this work could not be measured 
directly from the experimental data. For these we relied on Monte Carlo simulations of 
particle production and decay and of our detector's response to final state particles. The 
Monte Carlo program consisted of several parts: 

• We generated bb quark pairs according to the predictions of a next-to-leading order 
QCD calculation ]3"T| , |3"2] ] using the MRSDO parton distribution functions J53]. We 
required pr > 5 GeV/c for a 6-quark. We assumed the distribution in rapidity y to be 
flat in the range \y\ < 1.2. 

• We determined the b quark fragmentation into a B meson using the Peterson 
parametrization with the parameter e = 0.006 |S7]JS3 . 



For B c production we used the fragmentation model of Ref. [14 



• We used the CLEO B decay model |35|], for the decay of the B meson and its daughter 
particles. 

• We used full simulation of the CDF detector to calculate its response to the final state 
particles. 

The resulting Monte Carlo events were processed with the same programs used to reconstruct 
the data. The processes we studied with this program were: 

• B c -> J/ipe, 

• B c -> J/ipfi, 

• B —> J/ijj K, 

• Pairs of B mesons with B — > J/ipX accompanied by B — > e or \x either directly or 
through its daughters. 

These studies yielded ratios of the detection efficiencies e(B c — > J/ipe), s(B c — > J /if) //) and 
e(B — > J/tp K), the BB backgrounds described in Sec. |IV D| , and the K distributions used 
in Sec. |VT| 

In addition, we employed hybrid Monte Carlo calculations that replaced a real track 
in a. J/tp + track event by another particle to study punch-through, decay- in-flight, and 
photon-conversion backgrounds. These studies are described in Sec. |V|. 
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APPENDIX B: VALIDATION OF BACKGROUND ESTIMATES 



1. Semileptonic B Decay Sample 

We confirm our ability to determine accurately the various background rates to our 
observation of the B c meson by using identical methods to determine the background rate 
for a different process studied in a data sample independent of that which yielded the J/ip 
+ track distributions in Figs. |^ and 

In b hadron decays, leptons are produced either directly in the b —>■ c decay or in the 
sequential decay of the daughter charm hadron. Pairs of leptons thus arise from events in 
which there is a both a prompt and sequential semileptonic decay of a single B or from 
BB pairs. The leptons in the sequential decays are necessarily opposite charge and have a 
two-particle mass less than 5GeV/c 2 . Leptons from BB pairs may be of the same charge 
either because of mixing or where one lepton is direct and the second is sequential. The 
pair-mass, however, tends to be large and is typically greater than the B mass. Thus, low- 
mass, same-charge pairs of identified leptons in B events form a nearly pure background 
sample in which we can test our algorithms. 

Our overall strategy for obtaining such a sample was to select lepton pairs in which one 
lepton was responsible for the trigger and came from a displaced vertex. We required the 
other lepton also to originate in a displaced vertex in the same jet cone as the trigger lepton. 
This emphasized low mass pairs. 

Our inclusive, high-p^ lepton trigger provides a large sample of semileptonic b (and c) 
decays. However, even after strict identification cuts these events are contaminated by events 
in which the lepton is a misidentified hadron. Therefore, we need to identify the event as 
a B decay by other means. To do so, we take advantage of the long B lifetime. In central 
electron and muon events with lepton P T > 7.5 GeV/c, we reconstruct jets in the calorimeter 
using a cone algorithm |39| with a cone radius of R = y/rp + (j) 2 = 0.7. We require a jet of 



E T > 10 GeV and search for displaced decay vertices using charged particle tracks that lie 
inside the jet reconstruction cone. We define the impact parameter significance s = \do\/an 
where do is the impact parameter in the transverse plane with respect to the beamline, 
and ad is its measured uncertainty including the known transverse beam width. We require 
either that the lepton and two additional tracks in the cone satisfy s > 2.5 or that the 
lepton and one additional track satisfy s > 4.0. In all cases, we require that the displaced 
tracks originate from a common point and that the vertex be forward of the beamline with 
Lxy/cxy > 2.0, where a xy is the uncertainty on L xy . 

To estimate the purity of this sample, we make use of another property of semileptonic 
B decays. The lepton is typically the leading particle in the decay. Further, the lepton 
spectrum in the B rest-frame is well established [JRJ. In the candidate events, we find the 



distribution of the momentum of the lepton transverse to the jet direction Pr,rei and fit it to 
Monte Carlo templates for direct-6 and sequential decays, cc production, and false leptons 
from mismeasured prompt jets. We find a sample composition of approximately 85% bb, 
10% cc, and 5% false leptons. 

The tracks in the event, except for the trigger lepton, provide the parent sample to test 
the backgrounds to our soft-lepton identification. For each track that satisfies our electron 
or muon geometric requirements and comes from a displaced vertex in the same jet cone as 
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the trigger lepton, we find the mass of the trigger-lepton and candidate track combination. 
We weight the mass by the track's false lepton probability (as determined in section IV) 
and histogram the mass for same-charge and opposite-charge combinations. We compare 
this to combinations in which the candidate track satisfies our lepton identification criteria. 
Next-to-leading-order processes can contribute to the low-mass regions with leptons from 
different b hadrons. Therefore, to make an accurate comparison, we find the distribution of 
lepton-pair masses in BB Monte Carlo simulation subject to our trigger and identification 
criteria. We used the number of trigger leptons to normalize the BB Monte Carlo calculation 
to the experimental results. 

For various combinations of electrons and muons identified in the trigger and those 
identified in subsequent analysis (tagged) Fig. |27| shows the mass distributions of same- 
sign di-leptons. The points with uncertainties are the data, and the histograms represent 
the contributions from the same backgrounds relevant to the B c analysis. Table [V| lists 
the number of expected and observed di-lepton pairs for M« < 5GeV/c 2 . The calculated 
and observed same-sign di-lepton data are in reasonable agreement within the statistical 
uncertainties. This supports the validity of the background calculation in the B c analysis. 

We also removed the requirement that the second lepton come from a displaced vertex 
in the same jet cone as the trigger lepton and repeated the analysis with this larger sample. 
In this case, we normalized the BB contribution by requiring that the sum of same- and 
opposite-charge BB and false-lepton contributions in the high-mass (Mm > 5 GeV/c 2 ) region 
be equal to the total number of di-lepton events. The two normalization procedures agreed. 



2. Impact Parameter Significance 

We present additional evidence that the BB background, based on a Monte Carlo cal- 
culation, is indeed small. We re-analyzed the J/ip i data with a modified procedure which 
relaxed the requirements that the third track come from the same point as the J/ip decay 
vertex. 

• We performed a two-track mass and vertex constraint on J/ip — > fi + fi~ and required the 
good-fit probability to be greater than 1%. This departs from our standard procedure 
of requiring all three leptons to originiate at a common vertex. 

• With the third lepton, we calculated the J/i[)£ mass, px and ct* based on the J/ip 
vertex. 

• We required ct* to be greater than 60 fj,m. 

• We calculated the distance of closest approach of the third lepton track to the J/i[) 
vertex d and its uncertainty ad- We define the ratio djod as the impact parameter 
significance. 

Figure ^(a) shows the impact parameter significance for electrons with respect to a J/ip 
vertex for the data. Figure ^(b) shows the same quantity where the third lepton is a 
muon. Backgrounds from BB should extend to higher values of the impact parameter 
significance because the J/ip and the third lepton come from different vertices. B c events 
should populate the low impact parameter region because the J/ip and the third lepton 
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emerge from a common vertex. The figure shows that, when this region is included, most 
events have low impact parameters. Note that the events in Fig. 28 are a superset of our 
final data sample because of the relaxed vertex requirements. When we account for the effect 
of the relaxed requirements on these events, the level of events with high impact parameters 
is in good agreement with our predicted levels of BB backgrounds. 



37 



APPENDIX C: THE LIKELIHOOD FUNCTION 



For the likelihood analysis to test the null hypothesis and to estimate the size of the B c 
signal we used a normalized log-likelihood function. 

^-21n(|) (CI) 

where £ is the likelihood function, i. e. the product of all the probability distributions in the 
analysis, and Co is its value for a perfect fit. For purely Gaussian probability distributions, 
£ 2 is formally identical to the commonly used \ 2 • The advantage of £ 2 for a more general £ 
is that its properties are quantitatively similar to % 2 |] 

Below we define the input information and corresponding parameters along with the 
constraints and relationships among them, and we present the normalized log-likelihood 
function. Upper case letters represent input information, and lower case letters represent 
parameters of the fit. The superscript fi (e) refers to J/ipfi (J/ipe). We designate back- 
ground types by additional superscripts, fe and j \i for false leptons, ce for conversion 
electrons, Bfi (Be) for the BB contributions to muon (electron) backgrounds. We use Af 
(A?) to represent a function of the parameters corresponding to the fitted number of signal 
plus background events in the i th bin for the muon (electron) distribution. We use primes 
(N', n') for the number of events in the mass range 3.35-11.0 GeV/c 2 , and elsewhere we use 
unprimed numbers (N, n) for the subset in the range 4.0-6.0 GeV/c 2 . 

In order to propagate the uncertainties for various measured or calculated quantities, each 
item of input information has a corresponding parameter in the fit that we constrained to 
the measured value within its uncertainties. We include each such constraint as a Gaussian 
or Poisson factor in the likelihood function. For quantities with both Poisson statistical 
uncertainties and Gaussian systematic uncertainties, we adopted a Gaussian approximation 
of the Poisson uncertainty and added them in quadrature. 

It is important to understand that the only freely adjustable parameter in this fit is 
(a e + a M ), the total number of B c signal events. All other parameters are constrained within 
uncertainties by information independent of the B c candidate mass distribution. 

1. Definitions 

Data. For the histogram in Fig. |](b), we represent the number of candidate B c —>■ J/ip \i 
events in the i th bin as . These numbers contribute factors to C according to the Poisson 
probabilities 



6 As an example, if £ is a simple product of either Binomial or Poisson probabilities, it is easy to 
derive an expression for the inverse of the co- variance matrix for £ 2 in the same way one does for \ 2 ■ 
This yields the textbook uncertainties in the parameters. A Taylor expansion of the logarithmic 
terms in ^ 2 reveals that a one-standard-deviation change in a parameter from its best-fit value 
increases £ 2 by approximately one unit. 
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(A MU>: 



PvW) = ^t-* (G2) 



where the best estimate for the mean of Df is represented by Af, the function that sums 
the signal and background contributions calculated in the fit. Each term in the sum is a 
product of parameters defined below. In like manner, we symbolize the bin-by-bin numbers 
of candidate B c — > J /if} e events (Fig. |](b)) by Df and the functions representing their means 

B c Signal. The Monte Carlo simulation of B c production and decay, and response of 
our detector (App. [AJ) yielded mass distributions for J /if) \i and J/ipe and normalized each 
to unit area. Their values for the i th bin are represented by and Sf, respectively. We 
symbolize the total number of J/if>£ events by n' e and the fraction of these in the J/ipe 
channel by r e . For convenience, we express the numbers of events in the two decay channels 
as n ,fl = (1 — r £ )n' e and n' e = r £ n' e , and we emphasize that these are derived from the fitted 
parameters. The contributions to Af and A? are n'^S^ and n' e Sf, and their sum is shown in 
Fig. [l2|(a). The Monte Carlo simulation also determined the expected fraction of electron 
signal events, R e = 0.583 ± 0.043, which contributes a constraining Gaussian probability 
factor to C 

P(r e ) = — =J exp f- (r£ ~^n • (C3) 

False Muon Background. Fig. |](a) shows the mass distribution for the subset of 
J /if) + track events that satisfied the purely geometric criteria for third-track muons. The 
bin contents of this distribution are represented by This sample formed the parent 

distribution for calculating the false muon contributions from punch-through and decay- 
in-flight, and we combined these two sources of background into a single distribution. We 
calculated the bin-by-bin sums over the distributions in Figs. §(a) and (b) and normalized 
each distribution to unit area. To allow for a shape difference from the parent distribution, we 
calculated the bin-by-bin fraction F( M ±Ai</ M of the parent distribution. We also scaled these 
fractions so that the resulting products F/ M j/ M summed to 1.0. The quantities F/ M account 
for any shape difference between the parent J /if) + track distribution (j(^), and the false- 
lepton distribution (if" J* ), and they normalize the latter to unit area. N>f» = 11.4 ± 2.4 
is the total number of background events that satisfied all the muon identification criteria. 
In order to allow the fit to vary within the uncertainties in these measurements, we replaced 
them by parameters. The fitted parameters were constrained by Poisson contributions 
to C. There was a similar constraint for N'^. The fitted parameters were constrained 
through Gaussian factors in C The parameter, n'^, is constrained to N'^ L by Gaussian 
factor in C The contribution to Af from false muon backgrounds is n'^ f{ fl j{ fl ■ 

False Electron Background. Our background estimate for false electrons used another 
subset of the J /if) + track distribution that satisfied the purely geometric criteria for third- 
track electrons. This parent distribution is j/ e . (Below we discuss correlations between j/ M 
and j( e '.) The remaining input information and parameters for the false electron background 
are formally identical to those for the false muon background: F(\ f/ e , j/ e , N'f e = 4.2 ±0.4 
and n'f e . The contribution to Af from false electron backgrounds is n'f e ff e j{ e . 

Conversion Electron Background. We measured the number of identified conversion- 
electron background events to be N' ce = 2. We represent this by a parameter n' ce constrained 
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to N' ce by a Poisson factor in C. With the hybrid Monte Carlo calculation, we determined the 
ratio of residual (not identified) conversions to identified conversions to be R ce = 1.06 ±0.36. 
The corresponding fitted parameter is r ce constrained by a Gaussian factor in the likelihood 
function. We re-normalized the mass distribution for residual conversions in Fig. |(b) to unit 
area, represented by J 4 ce ± AJ 4 ce . The corresponding parameters are with constraining 
Gaussian probability factors in C. The contribution to A^ from residual conversions is 
n' ce r ce J_f e . 

BB Backgrounds. We used a Monte Carlo procedure to calculate independently the 
BB background contributions to J/ipfi and J/ipe. The shapes of these were found to 
be identical and a single parent distribution distribution, J B ± AJ B , normalized to unit 
area, was adopted for both. It is represented by the parameters jf that are constrained by 
Gaussian terms in the likelihood function. The Monte Carlo results for total numbers of 
events are: N' Be = 2.3 ±0.9 for J/ipe and N ,Bfl = 1.44 ±0.25 for J/ipfi. The corresponding 
parameters are n' Be and n' Bfl . The contribution to Af and A^ from BB background are, 
respectively, n' BtJ- j B and n' Be j B . 

Sums. We present here the two functions that, through their parameters, are adjusted 
for the best fit to the data distributions, and Df . 

Af = (1 - r>'^f + n"'"ff"jf" + n' B »j B (C4) 
\l = r £ n ,e S* + n' fe f/ e j{ e + n' ce r ce tf e + n' Be j B (C5) 



2. The Normalized Log-Likelihood Function 



It is easy to show that a likelihood function C, which is the product of factors of the 
form given in Eq. [C3| , leads to £ 2 = \ 2 through Eq. £JT]. F° r £ composed of Poisson factors 
like those in Eq. |C2| we find the corresponding factors in £ to be 



(D. 



--Df 



Dp. 



(C6) 



In the ratio P D ^(Xf)/P D ^(D^) the factorials cancel, and the contributions to the normalized 
log-likelihood function are quite simple. 



12 



-2 In 



2£ 



C 



Af 
D? 



(C7) 
(C8) 



where £' 2 is the first part £ 2 which we now write down in full. 



+ 



(Af - D?) - In 



Af 

J i 



+ 



(X-Bt)-Bt hi 



K 

Df 



J, 



fn 



+ 



fe' 



(C9a) 
(C9b) 
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Line |C9a| is the fit to the B c candidate distributions. Lines |C9b| and |C9c| constrain the 
parent distributions for the various backgrounds and the shape-dependent fractions for the 
false lepton distributions. Lines |C9d| and |C9e| constrain the normalizations for the five 
background distributions, the Monte Carlo calculation of the expected ratio of electron 
to muon B c events and the calculated ratio of residual to identified conversion-electron 
background events. 



3. Log-Likelihood for the Mass Analysis 

In the B c mass analysis, we performed an unbinned likelihood fit to the observed J /ip i 
mass distribution. The unbinned likelihood function for this analysis was the product of 
the probability distributions for the J/ipl mass for the B c signal and the backgrounds. The 
individual contributions to the probability distribution played a role similar to that defined 
for the bin fractions in App. |C 1| except that the bin index i was replaced by m^, the J/ipt 
mass for the i th event. Further, the signal distribution differs for each assumed B c mass 
M(B C ). 

• — > S^(rrii, Mb c ) and Sf — > S e (rrii, Mb c ) represent the normalized signal distribu- 
tions. 

• F^irrii) and F e (rrii) represent the normalized false /i and false e background distribu- 
tions. 

• J B {rrii) represents the distribution of the BB background obtained from Monte Carlo 
calculations. 

• J ce (rrii) represents the distribution for conversion and Dalitz decay electrons. 

The preliminary version of each of the above functions was as a smooth spline fit to the 
appropriate binned distribution. The fit was done prior to excluding events within 50 MeV 
of the B + mass to eliminate B — > J/ip K. The final version of the probability distribution 
was provided by a computer algorithm which, given a specific value for m«, returned the 
value of the spline function except when mj was within the excluded region for B — > J/ip K, 
in which case it returned zero. The areas of the final probability distributions were each 
normalized to unity over the range 3.35 to 11.0 GeV/c 2 . and N' e are the total numbers 
of /i events and e events. 

The normalized probabilities for the muon and electron distributions are X^/D^ and 
X e /D e , where 
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y( mi , M B J = (1 - r>'^K, M Bc ) + n'^F^rm) + n' B »J B { mi ) 
X e ( mj ,M Bc ) =r £ n' l S e {m,j,M Bc )+n' fe F e <- ^ 1 - ,Be lBe <- ^ 1 
D» = (1 - r £ )n /£ + n'^ + n'^ 



(nij) +n' Be J Be (m j ) +n' ce J c 



rrii 



(ClOa) 
(ClOb) 
(ClOc) 
(ClOd) 



The unbinned likelihood function contains the product of these probabilities for all the 
events. The parameters in the probability functions were adjusted for the best fit to the 
data. The likelihood function also contained constraints on the parameters determined 
independently of the candidate events. We define the log-likelihood function by 
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(C12b) 
(C12c) 

(C12d) 
(C12e) 



where C was chosen so that £4 = at £ = C m in- Line |C12a| is the fit to the B c candidate 



distributions. Line U12b is the constraint to the total numbers of J/ip fi and J/ip e events. 



Lines |C12c| and |C12d| constrain the ratio of e to /x signals and the number of background 
events for each background. 



4. Log-Likelihood for the Lifetime Analysis 



The unbinned likelihood function used to fit the B c lifetime was a product over the 42 
J/ipe and 29 J/ijj n candidates of the probability distribution for ct*. 

The normalized probabilities which combine both signal and background distributions 
in Xi = ct* for the J/ip /i and J/ip e are A^/D"^ and A e /D" e , where 



A^x u ct) = (1 - r £ )n" e ^ ig {xi, ct) + n" f ^ f » 



X ; 



n 



trBfi-pBfj,^ 



A e ( Xj , ct) = r £ n" l F e si { Xjl ct) + n" fe F fe 



JlBe T-Bei 



D" e 
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)J +n" ce F ce (x j 



+ n 



li Bp, 



r £ n" e + n" fe + n" Be + n" ce . 



(C13a) 
(C13b) 
(C13c) 
(C13d) 



The symbols are defined in Sec. |V11| . The A-functions, of course, depend on all the fitted pa- 
rameters, but we choose to emphasize the dependence on ct which is the only unconstrained 
parameter. 
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These probabilities are functions of the parameters given in Table |I| which describe 
the various backgrounds. For each background, the shape parameters were determined by 
a background fit that yielded the values in the table, which we represent by A J k , and a 
variance matrix V kl , where j is the background label and k and / label the three or four 
shape parameters. The lifetime fit varied a parameter, a 3 k corresponding to each of the A J k , 
and these were constrained by a contribution to the log-likelihood function 

x) = Ek - A i) ■ (vi)' 1 ■ w - A i) (° 14 ) 

k,l 

The number of events in each signal and backround contribution was subjected to a 
Gaussian or Poisson constraint as in the previous parts of this Appendix. 

We considered differences between the prompt background and that in the B-like region 
with ct* > 60 pm. Our dE/dx measurements indicated that for J/ip + track events, the 
pion fraction for the third tracks in the prompt region was f% = 74 ± 4% compared with 

= 56 ± 3.4% noted in Sec. [IV A 1| for the B-like region. These uncertainties are statistical 
only. In order to account for systematic uncertainties, we assumed po = f\j f% = 0.75±0.25. 
In the fit, we also allowed a variation in the relative probability ujq = 3.3 ±0.4 for pions and 
kaons to be mistakenly identified as electrons. The effect of this is to modify the values of 
some of f±, which become cumbersome algebraic functions of the fitting parameters p and 
u). For clarity in the equations, we omit these details. 

The log-likelihood functionf] that combines the unbinned fit to the ct* values for the 
candidate events and constraints on the parameters describing the probabilities is 

- 2 In C comb = -2 ln(£ e /^) (Cl5a) 
= -2^1nA e (xi) -2^1nA^(xi) (Cl5b) 

i i 

+2 [n" Bc + n" fe + n" ce r ce + n" Be + n" ffl + n" Bfl + \n(N" e \) + ln(N""\)] (C15c) 
+2{n ce - N" ce lnn ce + ln(A^ ,/ce !)) + ( T ~ Rce J (C15e) 



n fv — N"fi- 1 \ 2 / n B » — ]S[" B ^ \ 2 f n"f e — N"f e \ 2 ( n Be — N' 1 



jnBe \ 1 

+ { AN"fv ) + { AN" B ^ ) + ^ AN"f* ) + AN" Be ) ( ° 15f) 

, ( P~P0 \ 2 , f U-Wq \ 2 ,2,2,2 //"ii r \ 

+ {-Ap-) + {-A^) + X /e + X„ + X /M . (C15g) 

Note that terms N" e In D" e and N" 11 In D"^ do not appear because they cancel between the 
denominator of the log-probability sum (Line |C15b[ ) and the numerator of the Poisson con- 
straint on the numbers of J/ipe and J/ip p events (Line [CT5c]). Line |C15d| is the constraint 



7 The log-likelihood function used here has a minimum of -21n£ com6 = -430. This value depends 
on the fact that ct* was expressed in cm in the computer program, although we have used p,m in 
this report. Had the programs used pm, the value would have been higher by 2N" e ln(10 4 ) = 1308. 
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on the J/ipe fraction in the number of B c events. Line |C15ej contains the Poisson con- 
straint on the number of detected conversion electron background events and the Gaussian 
constraint on the ratio of undetected to detected background. Line |U151| contains Gaussian 
constraints on the numbers of other types of background events. Finally, Line |C15g| provides 
constraints on p, u, and the shape parameters for the background probability functions. 
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TABLES 



TABLE I. B c Signal and Background Summary: The Counting Experiment 

4.0 < M(J4>£) < 6.0 GeV/c 2 



J/ip e results J/ip \x results 



Misidentified leptons 
False Electrons 
Conversions 
Total False Muons 
Punch-through 
Decay-in-flight 
BB bkg. 


2.6 ±0.05 ±0.3 
1.2 ± 0.8 ± 0.4 

1.2 ±0.5 




6.4 ± 0.5 ± 1.3 
0.88 ±0.13 ±0.33 

5.5 ± 0.5 ± 1.3 
0.7 ±0.3 


Total Background a 


5.0 ± 1.1 




7.1 ± 1.5 


Events observed in data 


19 




12 


Net Signal 

Combined 


14.0 


18.9 


4.9 


Pcounting (Null) b 


2.1 x 10~ 5 




0.084 



a Upper limit on other backgrounds < 0.44. 

b The probability that the background can account for the data in the absence of of a signal is 
based on a convolution of Poisson uncertainties and Gaussian uncertainties of the backgrounds. 
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TABLE II. B c Signal and Background Summary: The Likelihood Analysis 



Input Constraint a 



(Results of Fit) a 




J/ip e results 




J/ip n results 


False Electrons 


N'f e = 4.2 ± 0.4 
(n'f e = 4.2 ± 0.4) 






Found Conversions 


N lce = 2 

(n /ce = 2.2 ± 1.4) 






Conversion ratio 


R ce = 1.06 ± 0.36 
(r ce = 1.08 ± 0.35) 






Unfound Conversions b 


2.1 ± 1.7 
(2.4 ± 1.7) 






False Muons 






N /fn = n 4 ± 2.4 
(n'^ = 9.2 ± 2.3) 


BB bkg. 


N' Be = 2.3 ± 0.9 




n ib^ = 1 44 ± 0>25 




(n' Be = 2.6 ± 0.9) 


[n lB ^ = 1.42 ± 0.25) 


Total Background 


8.6 ± 2.0 
(9.2 ±2.0) 




12.8 ±2.4 
(10.6 ±2.3) 


Total Signal 




(n« = 20.433) 




Electron Fraction 




R £ = 0.58 ± 0.04 
(r £ = 0.59 ± 0.04) 




e and fi Signal 


(n' e = 12.0^;») 




= 8.4+^) 


Signal + Background 


23 

(21.2 ±4.3) 




14 

(19.0 ±3.5) 


P(Null) c 




0.63 x 10~ 6 





a The numbers quoted here are for the mass range 3.35 < M(Jipt) < 11.0 GeV/c 2 . 
b Derived from other parameters. 

c Probability that background alone can fluctuate to produce an apparent signal of 20.4 events or 
more, based on simulation of statistical fluctuations. 
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TABLE III. Parameters for Background Distributions in ct* 



j fe fji ce Be Bfi 

N" j 13.2 ±1.3 12.6 ±2.8 a 1.5 ±1.1 0.79 ±0.34 

/| 0.199 ±0.004 0.36 ±0.01 0.45 ± 0.02 0.96 ± 0.01 0.98 ± 0.06 

fl 0.032 ±0.004 0.034 ±0.007 0.12 ±0.02 1 - /f e 1-/+^ 

(H 371 ± 15 445 ± 20 382 ± 27 371 ± 15 406 ± 16 

Ai (/im) 103 ±9 96 ± 16 138 ± 27 65 ± 15 48 ± 21 



a The number of conversion background events was calculated from identified conversions N" ce = 3 
and the ratio R ce = 1.06 ± 0.36. See App. [Tl] 
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TABLE IV. K(f.s.) 



a(ByBR(B^> TF) Derived From Various Experimental Searches 



Experiment 


final state (f.s.) 


K(f.B.) 


DELPHI a Ref. |18| 


J 111) 7T+ 

IT 


< (0.9 to 0.7) (90% C.L.) 




J/ipi + v 


< (0.5 to 0.4) (90% C.L.) 




J/tp 7r + 7T~7T + 


< 1.5 (90% C.L.) 


OPAL Ref. 


J /lb 7T+ 


< 0.6 (90% C.L.) 




J/4>al 


< 0.3 (90% C.L.) 






< 0.4 (90% C.L.) 


ALEPH Ref. @ 




< 0.2 (90% C.L.) 






< 0.3 (90% C.L.) 


CDF Ref. a 


J/lp 7T+ 


< (0.15 to 0.04) (95% C.L.) 


This Expt. 


J/ipl + v 


0.132l^ (stat.) ± 0.031 (syst.) ±8:8io (lifetime) 



a The ranges quoted for DELPHI and CDF Ref. [21| correspond to the assumed B c lifetime range 
0.4 to 1.4 ps. 
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TABLE V. Calculated and observed false leptons in the background validation 



Tagged Sample a 


e 




Observed e 


33 


37 


Expected background + BB 


43 ± 10 


38 ±7 


Observed fj, 


43 


63 


Expected background + BB 


41 ±4 


70 ±6 



a The numbers here are for events with dilepton mass < 5 GeV/c 2 . 
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FIGURES 




.85 2.9 2.95 3 3.05 3.1 3.15 3.2 3.25 3.3 3.35 



M(|iV") (GeV/c 2 ) 



FIG. 1. The distribution of masses. The data used for further analysis lie between 

3.047 and 3.147 GeV/c 2 and contain 196,000 ± 500 J/ip events above a background continuum of 
20, 000 ± 150 events under the J/tp peak. 
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FIG. 2. The distribution of masses of J j^K^ candidates. The solid curve represents a least 
squares fit to the data between 5.15 and 5.8 GeV/c 2 consisting of a Gaussian signal above a linear 
background. The area of the Gaussian contribution is 290 ±19 events. 
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FIG. 3. Histograms of the number of events vs. M(J/ip + track), (a) J/tpe candidates. 
For the 6530 events in this histogram, we assigned the electron mass to the third track and re- 
quired pt > 2.0 GeV/c. We applied the geometric criteria but not the particle identification 
criteria for electrons, (b) The 23-event subset of the distribution above that satisfies the electron 
identification criteria. Note that the bins in M(J/ip + track) are not uniform in width. The 
bin boundaries are indicated by tick marks at the top of the figures here and in subsequent mass 
histograms. The binning is discussed in the text. 
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FIG. 4. Histograms of the number of events vs. M(J/ip + track), (a) J /tp \i candidates. For the 
1055 events in this histogram, we assigned the muon mass to the third track and required pr > 3.0 
GeV/c. We applied the geometric criteria but not the particle identification criteria for muons. 
(b) The 14-event subset from (a) that satisfies the muon identification criteria. 
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FIG. 5. The transverse-momentum distribution for the J/tp£ system in B c candidates (line). 
It is compared with the normalized pt distribution for all backgrounds (dark shading) and with 
the pt distribution for B c — > J ' /ip tv events generated by Monte Carlo calculations (light shading). 
The latter is normalized to the fitted number of B c events determined in Sec. [V[ 
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FIG. 6. Mass histograms for backgrounds from hadrons misidentified as muons. (a) The sum 
of punch-through background contributions from 7r ± , K + and K~ . The dominant contribution to 
the punch-through background is from K + because of its lower interaction cross section, (b) The 
sum of decay-in-flight background contributions from -k^ and K^. (c) The contribution from 
BB background. These plots are normalized by their calculated contribution to the candidate 
distribution in Fig. |||. 
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FIG. 7. (a) and (b) show the p^-dependent probability for kaons and pions, respectively, to 
decay in flight and be misidentified as muons. The specific ionization dE/dx was used to determine 
the correct proportion of pions and kaons in the data. 
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FIG. 8. The probability of incorrectly identifying a hadron as an electron as a function of pr- 
Tracks from both the JET20 sample and minimum bias sample were used (App. |B|). (a) and (b) 
show the data for the isolation parameter / < 0.2 and / > 0.2, respectively. The probability 
averaged over the third-track momentum spectrum for the events in Fig. |3|a is (0.066 db 0.006)%. 
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M(JA|/+track) (GeV/c ) 



FIG. 9. The J/ipe mass distribution (a) for background events resulting from misidentified 
electrons, (b) For events in which the electron originated from a 7 conversion or Dalitz decay that 
was not identified as such, (c) For BB events in which the J/ip came from one parent and the 
electron from another. These plots are normalized by their expected contribution to the candidate 
distribution in Fig. |3|. 
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FIG. 10. The difference in dE/dx observed for the third track in J/ip + track events and 
that expected for an electron. For electrons, Q has a mean of Q e and a standard deviation of ctq. 
We scaled the difference to yield a distribution with zero mean and unit standard deviation for a 
pure sample of electrons, (a) The same events shown in Fig. ||(a), where we assigned the electron 
mass to thethird track, required pt > 2.0 GeV/c, and applied the geometric criteria, but not the 
particle identification criteria for electrons, (b) The subset of the distribution above that satisfy 
all the electron identification criteria except dE/dx. 
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FIG. 11. Conversion Background, (a) The momentum spectrum for the track from an elec- 
tron-positron pair that fulfilled the requirements for the third lepton and (b) the momentum 
spectrum for the other member of the pair. To obtain a larger data sample, we removed the 
requirements on ct* and P{x 2 ) for the vertex fit. In both (a) and (b) the data distributions are 
normalized to unit area. The Monte Carlo distributions are normalized to the data for px > 0.5 
GeV/c 2 . In the lowest bin in (b), the difference between data and Monte Carlo arises from the 
dropoff in track reconstruction efficiency for px < 0.5 GeV/c. This is the reason for the undetected 
conversion background. 
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FIG. 12. (a) A tri-lepton mass distribution for B c — > J/ijilv based on Monte Carlo calcula- 
tions. It is normalized to the fitted number of B c events. The distribution was generated under 
the assumption that the mass of the B c is 6.27 GeV/c 2 . There are negligible differences between 
the shapes for B c — > J/ipfi and B c — * J/tpe. Note that (93.0 ± 0.6)% of the area falls in the 
signal region 4.0-6.0 GeV/c 2 . (b) The normalized mass distribution for all backgrounds for both 
muon and electron channels, (c) The mass distribution for B c candidates in the data for both 
muon and electron channels. Note that each of these is a summary histogram, i.e. the sum of 
several individual histograms presented earlier. We emphasize that the fitting procedures use the 
full information from individual distributions rather than the sums. 
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FIG. 13. Histograms of the J/ij}£ mass that compare the signal and background contributions 
determined in the likelihood fit to (a) the data for J/ipe and (b) the data for J/ijjfj,. Note that 
the mass bins vary in width. 
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FIG. 14. Histogram of the J/ijj£ mass that compare the signal and background contributions 
determined in the likelihood fit to the combined data for J/ipe and J/tp^i. Note that the mass 



bins vary in width. The total B c contribution is 20.4^?'? events 
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FIG. 15. The variation of £ 2 = — 2 In (C/Cq) as a function of the number of B c mesons. For each 
fixed value of N(B C ) all other parameters were adjusted for the best fit. We find N(B C ) = 20. 4+55 
at the minimum. 
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FIG. 16. Each entry in this histogram is the result of a fit to a Monte Carlo pseudo-experiment 
that simulated the statistical properties of our data. The backgrounds were generated with the 
measured means and varied using Poisson or Gaussian statistics as appropriate. B c events were 
included with statistical fluctuations from the mean of 20.4 and bin-by-bin fluctuations. The 
resulting muon and electron events were fit as with the data. The values of £ 2 are histogrammed 
here and compared with the value found for the experimental data. 
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FIG. 17. Each entry in this histogram is the result of a fit to a Monte Carlo simulation of 
the statistical properties of this experiment. We generated the backgrounds randomly according 
to the measured means and varied using Poisson or Gaussian statistics as appropriate. The B c 
contribution was set to zero in generating the distribution. We then fit the resulting numbers 
of muon and electron events using the likelihood function. The fitting function included a B c 
contribution. The histogram above is a measure of the probability of finding a false B c contribution 
of size N(B C ) where none exists. Upward and downward fluctuations of the generated samples can 
require both positive and negative solutions for N(B C ). We chose to collect all negative solutions 
in the lowest bin in this figure where these events produce a prominent excess. The smooth curve 
represents a fit of a convenient extrapolation function (the sum of two Gaussians) to estimate the 
area beyond 20.4 events. 
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FIG. 18. Templates used to determine the quality of the fit to the mass spectrum for various 
assumed values of the B c mass. Of the 11 values used, four templates are shown here for the 
following values of M(B C ): (a) 5.52 GeV/c 2 , (b) 6.27 GeV/c 2 , (c) 7.27 GeV/c 2 , and (d) 7.52 
GeV/c 2 . In each case the histogram displays the binned results of the Monte Carlo calculation and 
the smooth curve is a spline fit to the histogram. 
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M(B c ) (GeV/c ) 

FIG. 19. (a) The relative log-likelihood function £^ from fits to the data for various val- 
ues of the assumed mass of the B c . Error bars on £^ represent its fluctuations with different 
Monte Carlo samples of B c events at the same mass. The parabolic curve is a fit to the plot- 
ted points with x 2 / n d.o.f. = 4.3/8. A horizontal line is drawn through the parabola's minimum 
which occurs at M{B C ) = 6.40 GeV/c 2 . Another line one unit above its minimum indicates the 
one-standard-deviation uncertainties of ±0.39 GeV/c 2 . (b) The fitted number of B c events vs. 
M(B C ). It is stable over the range of theoretical predictions for M(B C ), 6.1 to 6.5 GeV/c 2 . 
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FIG. 20. K = (M(B c )/pt(B c )/(M(J/i/} £)/p T (J/ip £)) distribution using Monte Carlo simula- 
tion (a) for the electron channel and (b) for the muon channel. 
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FIG. 21. Pseudo-proper decay length distributions for the background distributions, (a) 
J/tpe background from false electrons, (b) J/ifte background from conversion electrons, (c) J/ip fi 
background from false muons. (d) BB background. Its shape is the same for both J/tpe and 
J/ipfi. 
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ct (/xm) 

FIG. 22. Pseudo-proper decay length distributions for data with the fitted curve and 
contributions from backgrounds (a) for the electron channel and (b) the muon channel. 
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FIG. 23. Pseudo-proper decay length distribution for the combined J/ip /i and J/ipe data along 
with the fitted curve and contributions to it from signal and background 
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FIG. 24. The change in — 21n(£) from its minimum as a function of cr for the fit to the ct* 
distribution of B c candidates. 
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FIG. 25. Results from 500 pseudo-experiments to simulate the statistics in the B c lifetime anal- 
ysis: (a) — 21n£; (b) fitted lifetime; (c) solid (dashed) line for the negative (positive) uncertainty; 
(d) (cTjit — CTi n p U t) I 'a fa- In (d), the positive (negative) uncertanty was used when the fitted 
lifetime was smaller (larger) than the input lifetime. 
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FIG. 26. The point with 1-standard-deviation contour shows our measured value of the ratio 
a ■ BR(B+ J/ip l+X)/a ■ BR(B + J/ipK + ) plotted at the value we measure for the B c 
lifetime. The shaded region represents theoretical predictions and their uncertainty corridors for 
two different values of the semileptonic width T s j. based on Refs. [Q] and ||. The other numbers 
assumed in the theoretical predictions are discussed in the text. 
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FIG. 27. Same-charge di-lepton mass distributions for a trigger lepton and a tagged lepton. 
Both were required to come from a displaced vertex and be within the same jet cone. In (a) the 
tagged lepton is an electron, and in (b) the tagged lepton is a muon. In both cases, the data from 
trigger electrons and that from trigger muons are combined. The points with uncertainties are 
data, and the histograms show the predicted contributions from the various backgrounds relevant 
to the B c analysis. 
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(b) 4.0 < M(J/V m) < 6.0 
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FIG. 28. (a) Distribution of the impact parameter significance of the third track with respect 
to the J ftp vertex for J/ipe events, (b) The same distribution for the J/ipfx events. 
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